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The internal structure of the resonant A(1405) state is investigated based on meson-baryon 
coupled-channels chiral dynamics, by evaluating density distributions obtained from the form factors 
of the A(1405) state. The form factors are defined as an extension of the ordinary stable particles 
and are directly evaluated from the current-coupled meson-baryon scattering amplitude, paying at- 
tention to the charge conservation of the probe interactions. For the resonant A(1405) state we 
calculate the density distributions in two ways. One is on the pole position of the A(1405) in the 
complex energy plane, which evaluates the resonant A(1405) structure without contamination from 
nonresonant backgrounds, and another on the real energy axis around the A(1405) resonance energy, 
which may be achieved in experiments. Using several probe interactions and channel decomposition, 
we separate the various contributions to the internal structure of the A(1405). As a result, we find 
that the resonant A(1405) state is composed of widely spread K around A'', which gives dominant 
component inside the A(1405), with escaping ttE component. Furthermore, we consider KN bound 
state without decay channels, with which we can observe the internal structure of the bound state 
within real numbers. We also study the dependence of the form factors on the binding energy and 
meson mass. This verifies that the form factor defined through the current-coupled scattering am- 
plitude serves as a natural generalization of the form factor for the resonance state. The relation 
between the interaction strength and the meson mass shows that the physical kaon mass appears to 
be within the suitable range to form a molecular bound state with the nucleon through the chiral 
SU(3) interaction. 
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I. INTRODUCTION 



The A(1405) is a baryonic resonance state with spin- 
parity ~ 1/2^, isospin / = 0, and strangeness 
S = —1, and is located just below the threshold of anti- 
kaon (K) and nucleon (N). This resonance has been 
considered as a quasi-bound molecule state of the KN 
system [1-3] before the establishment of QCD. In sim- 
ple constituent quark models, the A(1405) is classified 
into the 70 dimensional representation of the spin-flavor 
SU(6) with excitation of one of the quarks to the p- 
state [4], but it was hard to explain the lighter mass of 
the A(1405) than the nucleon resonance iV(1535) in the 
same representation and the larger spin-orbit splitting 
between A(1405) and A(1520) than the nucleon part- 
ners. Recently the A(1405) was investigated with the 
unitarized coupled-channels method based on chiral dy- 
namics [5-10], in which the A(1405) is dynamically gen- 
erated in meson-baryon scattering without introducing 
explicit pole terms. In this method, one can successfully 
reproduce the cross sections of K~p to various channels 
together with the mass spectrum of the resonant A(1405) 
state below the KN threshold [5-8, 11-17]. 



composed of two resonance states [7, 10, 18].^ These two 
states have different coupling nature to the meson-baryon 
states [10, 18], and the A(1405) state which dominantly 
couples to the KN is located at 1420 MeV instead of the 
nominal 1405 MeV. These difference may be important 
for the KN system, as the A(1405) resonance position 
is measured from the KN threshold. The double-pole 
structure also suggests that the A(1405) resonance po- 
sition in the ttE invariant-mass spectra depends on the 
production mechanism of the A(1405). Some experimen- 
tal indications of this structure are found in Refs. [20-24]. 
Physical origin of the two poles is attributed to the at- 
tractive forces both in the KN and the ttS channels [25] . 

Recently, the importance of the meson-baryon dynam- 
ical component in the structure of the A(1405) has been 
discussed in the chiral unitary framework. It was revealed 
in an analysis of the phcnomenologically obtained scat- 
tering amplitude that the A(1405) is described predom- 
inantly by the meson-baryon dynamics [9]. The study 
of the A(1405) based on the Nc scaling in chiral uni- 
tary framework suggested the dominance of the non-qqq 
component in the A(1405) [26, 27]. The electromagnetic 
mean squared radii were shown to be much larger than 



One of the important consequences on the structure 
in the coupled-channels approach is that the A(1405) is 



^ The two-state nature of A(1405) was pointed out first in Ref. [19] 
in a difl'crent model. 
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that of ground state baryons [28].^ Relationship between 
the couphngs to the channels and the wave function of 
the A(1405) was clarified in Ref. [30]. The spatial struc- 
ture of the A(1405) was also discussed based on the KN 
molecular picture with the chiral KN potential [31, 32]. 
The molecular-like structure of the A(1405) suggests fur- 
ther few-body nuclear systems with kaon [33], such as 
IvNN [31-39], KKN [40-42], and KKN [43]. Although 
the dominant component of the A(1405) was found to 
be the meson-baryon molecule structure in these studies, 
the internal structure of the A(1405) resonance in terms 
of the hadronic constituents was not directly extracted 
from the meson-baryon scattering amplitude. 

If the A(1405) is dominated by the meson-baryon 
quasi-bound molecule, it may have a spatially larger size, 
which will lead to the softer form factor than that of typi- 
cal baryons dominated by the genuine quark component. 
The soft form factor could result in some experimental 
consequences, for instance in the nontrivial energy de- 
pendence of the A(1405) photoproduction cross section 
recently reported in Ref. [44]. The information of the 
size of hadrons is also important for the production yield 
in the heavy ion collisions estimated by the coalescence 
model [45]. One of the standard methods to study the 
structure of a particle is to use the vector currents, such 
as the photon, as probes. The properties of the A(1405) 
in relation to electromagnetic dynamics have been inves- 
tigated in the chiral unitary approach in Refs. [18, 46- 
51]. In the molecular picture of the A(1405), the pho- 
ton coupling to the A(1405) is introduced through the 
photon couplings to its constituent meson and baryon. 
Among others, the form factor of the A(1405) will tell 
us about the intuitive "size" of the particle. Therefore, 
for both theoretical and experimental understandings of 
the A(1405), it is interesting to investigate the internal 
structure of the A(1405) through the evaluation of its 
form factor. 

In this paper we study the form factor of the A(1405) 
in chiral unitary approach, in which the A(1405) is de- 
scribed as the dynamically generated resonance state. 
We emphasize that this is the first study to extract the 
form factor of the A(1405) directly from the scattering 
amplitude involving the resonance state as a response to 
the external probe current in a microscopic way. We have 
already reported the result of the electromagnetic mean 
squared radii of the A(1405) in Ref. [28]. Here we extend 
this scheme to evaluate the momentum dependence of the 
form factors, using the chiral effective theory of the octet 
mesons and baryons with the external currents for the 
elementary interaction between the constituent hadrons 
and the currents. From the form factors the density dis- 
tributions can be obtained through the Fourier transfor- 



^ The electromagnetic mean squared radii of the A(1405) was eval- 
uated also in [29] within the framework of the bound-state soliton 
model. 



mation. This allows us to visualize the spatial distribu- 
tion of the A(1405) in coordinate space. We also include 
the finite size effect of the constituent hadrons. 

There is one subtlety in the formulation because the 
A(1405) is a hadron resonance and decays via strong 
interaction. In this case, the form factors (or coupling 
constants in general) are obtained as complex numbers, 
whose interpretation is not as straightforward as in the 
case of the stable particle. Such a difficulty about the 
complex form factors generally shows up for the unstable 
states, as seen in the A(1232) electromagnetic form fac- 
tors studied in Refs. [52, 53] as complex numbers. Here 
we carefully define the form factor of the resonance par- 
ticle as an extension of the ordinary stable particles, and 
explain the methods to extract the form factors out of the 
photon-coupled meson-baryon scattering amplitude in a 
gauge invariant way. In order to achieve the intuitive 
understanding of the size of the resonance, we evaluate 
the form factor with several different model set-ups and 
draw a conclusion through the comparison of the results. 

This paper is organized as follows. In Sec. II we ex- 
plain our scheme to calculate the form factors and the 
density distributions of the A(1405) based on the chiral 
dynamics, taking account of how to extract the informa- 
tion of the "size" of the resonance state. The coupling 
of the A(1405) to the external current is obtained by 
the couplings of its constituent meson and baryon. In 
Sec. HI we show our numerical results of the calculations 
for the resonant A(1405) state, on the resonance pole po- 
sition in the complex energy plane and on the real axis. 
To obtain the intuitive understanding of the A(1405) as 
a quasi-bound state, we consider the KN bound state 
without decay width and study its structure by varying 
the parameters such as the binding energy and the meson 
mass in Sec. IV. Section V is devoted to the conclusion 
of this study. 



II. FORMULATION 

A. How to probe structure of resonance state 

First of all, we discuss the method to probe the prop- 
erties of a short-living resonance state. Let us consider 
an observable represented by a Hermitian operator O. 
One of the simplest and model-independent ways to de- 
fine resonance properties is given by the evaluation of the 
matrix element of the operator O for the resonance state 
vectors 

(Zr(p')^"^PI^r(^)^+^ (1) 

This is an extension of the matrix element of the stable 
particle such as nucleon, (iV(-) 1017V(+)). The labels "-h" 
and "— " represent "in" and "out" states, respectively. 
Here the kct 1Zr(P)(+)) and bra (Zr(P)(-)1 vectors 
represent the resonance state in the generalized Hilbert 
space, whose eigenvalues of the energy-momentum opcr- 
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ator arc complex as^, 

p^|Zr(p)(+)) = (£;r„ p)IZr(p)W) 

(ZR(P)(-'|P^ = (i?R, P) (Zr(P)(-)| 



(2) 
(3) 

(4) 



Here Mr and Fr are the mass and decay width of the 
resonance state, respectively. It is noted that this def- 
inition leads to complex values of the matrix elements 
such as the radius and charge distribution for the reso- 
nance state, because the Hermitian operators are allowed 
to have complex eigenvalues for the resonance state in 
the generalized Hilbcrt space. When the decay width 
of the resonance state is sufficiently small, the imagi- 
nary part of the eigenvalues of the Hermitian operators 
is also small and the matrix elements can be interpreted 
as physical osbervables in the ordinary sence. In this pa- 
per, we mainly discuss the real part of the eigenvalues 
of the Hermitian operators as "observables" , assuming 
that the significance of the imaginary part is sufficiently 
small. 

In order to scan the internal structure of the resonance 
with different probes, we evaluate the matrix elements 
of various conserved vector currents, the electromag- 
netic transition current Jg^' baryonic current Jg and 
strangeness current Jg , stemming from the U(l) symme- 
try of the underlying theory. Electromagnetic current is 
obtained from the gauged U(l) symmetry with appropri- 
ate charges of quarks, baryonic current from the simul- 
taneous U(l) phase transformation for all the quarks, 
and the strangeness current from the phase rotation for 
strange quark only. The other currents may also be used 
as the probes of the structure, as long as the currents are 
well defined. 

The structure of the resonance state is reflected in the 
form factors, as in the same way with stable particles such 
as the nucleon. For resonance states with baryon number 
1 and spin 1/2, we define the form factors F^^{Q'^) as the 
matrix elements of the probe current J^{x), 

{Z^{P'f~'^\J>^{x)\Z^{P)'^+^) = e'^(^'-^)"--F^(Q2), 

(5) 

where Q2 ^ -{P' - PY{P' - P) ^. Especially in the Breit 
frame the matrix element can be written as [18]: 



(ZR(q/2)(-)|J''(x)|ZR(-g/2)(+)) 



^ ( -^time('9 )j -Pspacc(Q ) 



2, »<T X q 

2Afn 



Xz^, (6) 



with the form factors of the time component Ftinic(Q^) 
and the space one i^space(Q^)j the momentum transfer 



^ For the thorough description of the eigenvector for the resonance 
state, see, e.g., Ref. [54]. See also Ref. [55]. 



= {0,q), = , the spinor for the resonance state 
X^Ij^, and the Pauli matrices cr" (a — 1, 2, 3) for the 
spin space. For the electromagnetic current, the time 
component form factor is identified as the electric form 
factor Fe{Q^), whereas the space component form fac- 
tor as the magnetic form factor Fm(Q^). The magnetic 
form factor i^M(Q^) is normalized as the nuclear magne- 
ton /xn = e/{2Mp) with the proton charge e and mass 
Alp. We define the baryonic form factor Fb(Q^) and the 
strangeness form factor Fs{Q'^) as the time component 
of the matrix elements of the baryonic and strangeness 
currents, respectively. 

The form factors F{Q^) contain the information of the 
structure of the system in response to the probe interac- 
tion, being related to the "classical density" of the system 
(such as charge density) via the Fourier transformation 
as. 



p{r) 



-F{Q% 



(7) 



(27r)3 E{Q^) 

where r is the radial coordinate from t he center of mass 
of the system, and M and E{Q^) = ^/M^+Q^ are the 
mass and the energy of the state, respectively. Taking 
the nonrelativistic reduction {M/E 1), we have. 



p{r) 



(27r)3 



We can also obtain the mean squared radius: 



dF 



drrp{r)=-Q—^ 



(8) 



(9) 



For the magnetic mean squared radius, we take conven- 
tional normalization by dividing the right-hand side of 
Eq. (9) by the factor Fm(0). The mean squared radius is 
well defined when the density p{r) sufficiently falls off at 
large r and the integration of r"^ p{r) converges, which is 
the case on the poles for the bound states as well as the 
resonance states. 



B. Matrix elements of resonance state in scattering 
amplitudes 

Here we explain the method to extract the matrix el- 
ements of the resonance state from the scattering ampli- 
tudes of dynamical approaches [18, 28]. Once the matrix 
elements are extracted, following the prescription given 
in the preceding subsection, we can evaluate the form 
factors [from Eq. (6)] and density distributions [from 
Eq. (8)]. 

Let us here consider s-wave meson-baryon scattering 
[MB M'B') amphtude T[y/s) with the total energy 
i/s, in which the excited baryon appears as a resonance 
state. Performing the analytic continuation of the energy 
variable to the complex plane -^/s — >■ z, we obtain the 
scattering amplitude in the complex energy plane. The 
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(a) r,, 



(b) Ti;, 



FIG. 1: (a) Feynman diagram for meson-baryon scattering 
amplitudes close to the pole of the excited baryon, (b) Fey- 
namn diagram for meson-baryon scattering amplitudes with 
the probe current attached to the excited baryon. The dou- 
ble, dashed, solid and wiggly lines correspond to the excited 
baryon, the ground state meson, the ground state baryon and 
the probe current, respectively. 



^ig,{2TT)H\P^' ~ P'^'), (14) 

where |$i(P')^"'') is the out (in) s-wave 

meson-baryon scattering state with channel i (j) and 
momentum P'^ (P^)- The total energies are given by 
z" = y/P"''P^ and z = ^/WP'^. In this way, the residue 
of the pole in the scattering amplitude is related to the 
inner product of the resonance state with the scattering 
state, representing the coupling strength. It is important 
that 



FIG. 2: Examples of the Feynman diagrams for the scattering 
amplitudes TJf in which the probe current does not couple to 
the excited baryon. The double, dashed, solid and wiggly lines 
correspond to the excited baryon, the ground state meson, the 
ground state baryon and the probe current, respectively. 



resonance state is expressed as a pole of the scattering 
amplitude in the second Riemann sheet. We define the 
"pole contribution" of the amplitude as 



-iT,j{z)\ 



pole 



«5» 



Z - ^R, 



(10) 



where the pole position is expressed as the resonance 
mass and width as = Mr — iTn/'i and the indices i 
and j represent the final and initial meson-baryon chan- 
nels, respectively. The residue giQj at the pole posi- 
tion can be interpreted as the product of the coupling 
strengths to the meson-baryon channels i and j. On 
the real energy axis, the full scattering amplitude can 
be written as 



-iT,j{y/s) = igi 



(11) 



The contributions which are not represented by the pole 
term is summarized as a nonresonant background term 
Tf^ (see also [56]). 

Next we relate the meson-baryon scattering amplitude 
with the matrix element of the resonance state vectors 
introduced in the preceding subsection. The pole contri- 
bution in the scattering matrix is written as 

pole 



d^P" 



^($,(p')<~'I^r(^")*+') 



x(Zr(p")(-)|$,(p)^+)) 

-i(27r)^(54(P"-P'^)- 



(12) 
(13) 



P 



2r 



_rf*p_ 



(+)) 



(-)l 



(15) 



in Eq. (12) can be understood as the "projection oper- 
ator" to the resonance state of the mass Zr , which is a 
generalization of that of the usual stable particle. This 
resonance component of the amplitude, Tij(z)|poic, is rep- 
resented diagrammatically in Fig. 1(a). 

In the same way, we also consider the scattering am- 
plitude T^,j for the MBj* M'B' process, where 7* 
stands for a probe such as photon for the electromagnetic 
current. The Fourier component of the matrix element 
of the probe current in the s-wave meson-baryon state, 
5";^^^ , is given as: 

S^.^iP', P; Q2) 

= / d^xe-'i"'="{^,{P')^-^\iJ^{x)\'Pj{P)'^+^), (16) 



with incoming and outgoing momenta of the meson- 
baryon system P^ and P'^, respectively, and = 
—q^q^. The resonance contribution to the matrix el- 
ement is obtained by inserting the "projection opera- 
tor" (15) for the resonance state on both sides of the 
current operator: 



(27r)4 J (27r) 



X {'i>,{P'Y-^\ZR{P"'Y+^) 



z'" - Z^ 



X (Zr,(P"')(-'|*J^|Zr(P")<^ 

X (ZR(p")(-'|a>,(p)(+)) 



z" - Zt, 



z - Zr z — Zr 

(17) 

where we have used Eqs. (5) and (14). Hence, the pole 
contribution to the scattering amplitude of A/P7* — >• 
M'B' process is 

T^,,(z', z; Q2)| ^ _g^^_F{Q^)^—g^, (18) 

z' - Zr z - Zr 
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which is represented diagrammatically in Fig. 1(b) (here 
and following in this subsection we omit the superscript fi 
in Ty and F). It is important that this term contains the 
matrix element of the probe current with the resonance 
state as the residue of double pole a.t z = z' = Zji. 

The rest contributions of the amplitude T^ij has less 
singularity at the resonance position Zr. We show sev- 
eral examples of these contributions diagrammatically in 
Fig. 2, in which the current does not couple to the inter- 
mediate resonance propagator. 

Combining Eqs. (10) and (18), the matrix elements in 
the Breit frame can be evaluated as the residue of the 
single pole at z = z' = of —Tjij/Tij in the complex 
energy plane: 



F{Q'^)\Brcit — 



(z' - Zu)Ty,,{z', z- Q2) 



T^,{z) 



(19) 



As one can see from the above derivation, this scheme 
indicates that only the amplitude Tjj (T^ij) with the 
single- (double-) pole terms contributes to the right-hand 
side, whereas both nonresonant background term T^'^ 
in Eq. (11) for Tij{z) and the less singular terms for 
Tjij{z\ z; Q^) automatically drop and have no effect to 
the form factors for the resonance state. As confirmed in 
Rcf. [28], the form factor obtained in Eq. (19) is a gauge 
invariant quantity. 

Equation (19) should be evaluated in the complex en- 
ergy plane for the resonance state. It is an interesting 
issue how the matrix elements of the currents for the res- 
onance can be obtained on the real axis, in which ampli- 
tudes are in principle observed experimentally. In order 
to keep closer connection to experimental measurements, 
we discuss a method to evaluate the form factor of the 
resonance state on the real axis, as developed in Ref. [18]. 

As we have seen above, the matrix element of the cur- 
rent for the resonance state is expressed as the residue 
of the amplitude T^{z' , z\ Q^) at the resonance pole po- 
sition (z' = 2 = Zji). If we take the Breit frame of the 
resonance (z' = z). the form factor F{Q^) of the reso- 
nance can be written on the real axis as. 



T 



9i 



-^L^+T\%^{^s) (20) 



where T}^°^^ represents less singular terms than the 
double-pole contribution given in the first term. It is im- 
portant to note that, if the decay width of the resonance 
state is small and its pole position is close to the real 
energy axis, the double-pole amplitudes will give domi- 
nant contribution to the MB^* — >■ M'B' process than 
the other contributions in the resonance energy region. 

Motivated by Eq. (20), we define an effective form fac- 
tor on the real axis as 



F^«'(Q^ 



(21) 



As discussed in Appendix A, the effective form factor (21) 
is correctly normalized so as to give the corresponding 
charge of the system at = independently of \/s, if 
we consider all of the appropriate diagrams giving both 
the double-pole and less singular terms. 

Let us study the relation between the effective form 
factor (21) on the real axis and the form factor at the 
pole position (19). Around the resonance energy region, 
where the pole contribution dominates the amplitude, 
the derivative of the MB scattering amplitude (11) with 
respect to yfs is written as 



(a/s - Z-R 



BG 



(22) 



Assuming that the background contribution is smoothly 
changing with respect to the energy ^/s around the res- 
onance energy and, thus, neglecting the second term of 
Eq. (22), we can write the effective form factor approxi- 
mately as 

F^^'iQ'; V^) « F{Q') Tl^W^ ^'^"^ . (23) 

In this way, the effective form factor is related to the res- 
onance form factor F{Q^). The deviation mainly comes 
from the less singular terms in T^ij which contribute to 
the effective form factor F'^^. Nevertheless, as seen in 
Eq. (23), if we take the energy y/s close to the resonance 
mass and choose the channel which strongly couples to 
the resonance, we can reduce the contamination from 
the less singular terms. In addition F°^{Q'^) will coin- 
cide with F{Q^) evaluated by Eq. (19) when we take 
^/s by analytic continuation, since in this case 

both the less singular contribution T^°*5>^ in Eq. (23) and 
the nonresonant background contribution in Eq. (22) au- 
tomatically drop at — > Zr. 



C. A(1405) in chiral dynamics 

In this subsection, we briefly review our formula- 
tion of the A(1405) resonance generated dynamically in 
the MB — > M'B' process using chiral unitary model 
(ChUM) developed in Refs. [5-8]. It turns out that ex- 
plicit pole terms are not necessary in the elementary in- 
teraction for the description of the A(1405) in ChUM [9]. 

The starting point of our formulation is the important 
fact that chiral symmetry and its spontaneous breakdown 
in QCD constrain the form of low-energy interactions in- 
cluding Nambu-Goldstone bosons. This is systematically 
expressed in chiral perturbation theory [57-59] , in which 
the effective Lagrangian is sorted out according to chiral 
expansion. From the lowest order meson-baryon chiral 
Lagrangian, the tree-level s-wave meson-baryon interac- 
tion known as the Weinberg- Tomozawa term can be ob- 
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tained as, 



where, 



M, - M, 



(24) 
(25) 



with the incoming and outgoing meson momenta, and 
k'^, the meson decay constant /, the Clcbsch-Gordan 
coefficient Cij fixed by the SU(3) group structure of the 
interaction, and masses of outgoing and incoming baryon. 
Mi and Mj, respectively. The last form in Eq. (25) is 
obtained by applying the nonrelativistic reduction for the 
baryons. The explicit value of the coefficient Cij for the 
KN scattering is given in Ref. [6]. 

Only this lowest-order Weinberg- Tomozawa interac- 
tion, however, is not sufficient for the description of the 
scattering amplitude, especially for three flavors. The 
main reason is the existence of the baryonic resonance 
state A(1405) in / = and S ~ —1 channel, just be- 
low the KN threshold, which spoils any perturbative 
expansion around the threshold. Therefore, in order to 
reproduce the KN scattering amplitude, some nonper- 
turbative and coupled-channels treatment for the KN 
scattering is needed. 

One valuable approach to take into account the non- 
perturbative effect is to formulate a scattering amplitude 
fulfilling exact unitarity with the N/D method [60]. Fol- 
lowing Ref. [7], with assumption that the intermediate 
states of the scattering are composed of only one octet 
meson and one octet baryon and with neglect of the left- 
hand cuts, one can write the inverse of the general form 
of the scattering amplitudes fulfilling unitarity as, 



(26) 



s- Sq 
27r 



Pt{s) 



ds'- 

+ (s' - s)(s' - So) ' 
(27) 

(28) 



Here sq denotes the subtraction point, sf is the threshold 
value of s of the channel i and qi center-of-mass momen- 
tum in the channel i 



{s - Mf + mf )2 - Asm: 
4s 



(29) 



where m; represents the meson mass in the channel i. 

The function Gi takes the same form as, except for an 
infinite constant, the ordinary meson-baryon loop inte- 
gral. 



2Mi 



(27r)4 ql - mf + ie{P- qif - Mf 



(30) 

where e is an infinitesimal real constant to specify the 
boundary condition. With dimensional regularization 
scheme, which keeps analytic properties of the loop func- 
tion, the finite part of the loop integral Gi can be written 
as. 



G,(yi) = 



2M, 
167r2 



ai(Mr 



In 



Ml 



' Mf + s ^ 
In — 

2s M? 



2(?,Vi) + ln(s + Mf -m^ + 2q,yjG) 



ln(-s + Mf 



2q,yG) - ln(-s - Mf + : 



2q, 



(31) 



with the regularization scale /ireg and the subtraction 
constant ai = ~{i7T)'^ai/{2Mi). 

For the Tij in Eq. (26), we adopt the matching scheme 
developed in Ref. [7], which tells that the identification 
of 7ij with the tree-level amplitudes by chiral perturba- 
tion theory is valid up to 0{p^). Hence in our approach 
we identify the Weinberg- Tomozawa interaction Vij in 
Eq. (25) as Tij. In this way, the T-matrix is written in 



matrix form as, 

r(V^) - - G]-\ (32) 

This amplitude satisfy the following equation 

k k 

(33) 
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which corresponds to the Bethc-Salpeter equation in al- 
gebraic form. We hence refer to this as the BS ampli- 
tude. In the present approach we have 10 meson-baryon 
channels, i,j = K'p, K°n, 7r°A, 7r°E°, 77A , rjT,^, tt+S", 
7r-S+, K+E-, and K°E^ . 

Now we fix the parameters in our approach; we have 
the masses of the ground state mesons and baryons, the 
meson decay constant /, and the subtraction constant 
ai(/^rog) with the regularization scale ^reg- We use the 
isospin-averaged masses for the mesons and baryons: 



TABLE I: Coupling strengths of the A(1405) resonances, Zi = 
1391 - 66i MeV (upper) and Z2 = 1426 - 17j MeV (lower), 
to the meson-baryon state obtained as the residues of the 
meson-baryon scattering in the particle basis. 



Channel 


coupling [Z-i] 


KN {K-p, K°n) 


-0.86 + 1.26i 


ttE (7r°E°, TT+E", 7r'E+) 


-1.42-h0.88i 


TjA 


-0.01 -h0.79i 


KE {K°E°, K+E-) 


-0.33 -H0.30i 



Mn = A/p = A/„ 
Ma 

Afs = M^o = = M^- 



= 138.04 MeV, 
= 495.67 MeV, 
= 547.45 MeV, 
= 938.92 MeV, 
=1115.68 MeV, 
=1193.12 MeV, 
=1318.11 MeV. 



For the meson decay constant /, we choose an averaged 
value / = 1.123/^ with = 93.0 MeV, which is one of 
the typical values used in ChUM. Finally we choose the 
subtraction constant a,; so as to reproduce the threshold 
properties of K^p observed in stopped capture on 
hydrogen [61, 62], as done in Ref. [5], 



Channel 



coupling {Z2) 



{-p, K°n) 



KN {K-p, 
ttE (7r°E°, TT+E", 

TjA 

KE {K°E\ K+E- 



-E+) 



1.84 + 0.67i 
0.26-H0.85i 
1.44-h0.21i 
0.09-H0.24i 



Zi strongly couples to the ttE state, whereas the higher 
pole Z2 dominantly couples to the KN state. Among the 
poles, the higher one is considered to be originated from 
the KN bound state [25]. Since we are interested in the 
structure of the A(1405) in the KN bound state picture, 
we mainly consider the contribution from the higher pole 
Z2 and regard the lower pole Zi as the background. 



a^jv = -l-84, a.s = -2.00, a,A = -1.83, ^^^^ 
flrjA = —2.25, ariY. ~ —2.38, a^H = —2.67, 

with the regularization scale /iicg = 630 MeV [11]. 

In the present model, the A(1405) is dynamically gen- 
erated in the obtained BS scattering amplitude. With 
above parameters, the model well reproduces the A(1405) 
mass spectrum below the KN threshold [23]. Further- 
more, in the present model the A(1405) is expressed 
by two poles of the scattering amplitude in complex 
energy plane, as {Zi = 1391 — 66z MeV) and {Z2 = 
1426 - m MeV) [10]. The study of the renormalization 
condition reveals that the A(1405) resonance in this ap- 
proach is dominated by the meson-baryon molecule state 
and the effect from the possible seed of the resonance 
in the subtraction constant is found to be small [9]. It 
has been reported in Refs. [17, 25] that the position of 
the lower pole Zi is largely dependent on details of the 
model parameters, whereas that of the higher pole Z2 
shows little dependence. 

The residues of the meson-baryon scattering amplitude 
Tij(^/s) at the pole position express coupling strengths 
of the resonance to the meson-baryon channels as seen 
in Eq. (10). The coupling strengths of the two poles Zi 
and Z2 obtained in the present formulation are listed in 
Table I."* From this Table, one can see that the lower pole 



D. Electromagnetic interactions in chiral dynamics 

As discussed in Sec. II B, photon-coupling to the reso- 
nance state is obtained as the residue of the second-rank 
pole in the MBj* M'B' amplitude. This process can 
be calculated by making a photon couple to the scat- 
tering process MB M'B' which has been obtained 
in the previous subsection to describe the A(1405). For 
the calculation of the MBj* — >■ M'B' process, we take 
a picture that the photon couples to the resonance state 
only through its constituent mesons and baryons [18, 28], 
which may be valid in view of the dominance of the 
meson-baryon component in the A(1405) [9]. Thus, we 
need only electromagnetic interactions of the constituent 
mesons and baryons in our approach. 

Due to the requirement of the gauge invariance, the 
elementary couplings of the photon to the mesons and 
baryons should be given by gauging the chiral effective 
Lagrangian in a consistent way with the description of 
the A(1405). The photon couplings to the meson and 
baryon appearing in the BS amplitude are derived in the 
minimal coupling scheme. In addition, the anomalous 
magnetic couplings in the BB'^ and MBM' B'j vertices 
are given by the chiral perturbation theory as done in 
Ref. [18]. The sum of the above two contributions deter- 



* The coupling strengths are shown in the particle basis. To 
compare with the coupling strengths from the amplitude in the 



isospin basis as in Ref. [10], a factor y'N; should be multiplied 
with the isospin multiplicity Nj = 2 for KN and KE,, and Nj = 3 
for vrE. 
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mines the elementary electromagnetic couplings V^, for 
MM'7, V^. for BB'-f, and Tf^- for MBM'B'-f, respec- 
tively. 

Now let us consider the minimal coupling scheme. The 
photon coupling to the meson, Vp^ . . is given by 



'tV^S^^ k') ^ iQmXk + k'Y , 



(35) 



with the meson charge Qyii , the incoming and outgoing 
meson momenta and k'^. The minimal coupling of 
the photon to the baryon is given by 



(N), 



\P, P') 



2M, 



P 



2M, 



icr X q 

2M„ 



, (36) 



with the baryon charge Qb^ and the incoming and out- 
going baryon momenta p^^ and p''^. These two cou- 
plings (35) and (36) are appropriate with the propagators 
in the loop function (31). For the MBM' B'j coupling, 
we use the following vertex, which is required by the 
Ward-Takahashi identity at tree-level with the Weinberg- 
Tomozawa interaction (25): 



iV 



(N),M 



pin 



with the incoming and outgoing meson-baryon total mo- 
menta and P'^, respectively, and = Qvii + Qb^- 
We note that this MBM'B'-f coupling does not con- 
tain the magnetic part. Actually for the electromag- 
netic properties of the neutral excited baryon this term 
does not contribute due to Qem = Qm + Qb = 0. This 
is a different point compared with Ref. [18], where the 
normal magnetic part coming from the original form 
of the Weinberg- Tomozawa term (24), proportional to 
QMi + Qmj , was introduced. 

For the anomalous BB'j and MBM' B'^ couplings 
which are gauge invariant by themselves, we use the inter- 
action Lagrangian appearing in the chiral perturbation 
theory [63]: 



AM, 



■b^Tr{B[S>^,S''][F+B] 



^6°Tr(S[5^5^{F+,i?}) 



with 



4M 



F^^ = -e ( u^QFf^yU + uQF^yW 



(38) 



(39) 



vertices (a = 1, 2, 3): 



^' ^ 2Mn 



(40) 
(41) 



where we have made nonrelativistic reduction. Here the 
anomalous magnetic moment for the baryon Ki and the 
matrix Aij are given by. 



K, = h^d, + hlU 



A,. 



2/2 



(42) 



(43) 



where the coefficients d^, /i, Xy , and Yij are fixed by 
the fiavor SU(3) symmetry and their explicit values are 
found in Ref. [18]. The values of the coefficients and 
&g are determined to be = 2.40 and &g = 0.82 so 
as to reproduce the observed anomalous magnetic mo- 
ments of the baryons^. In the calculation, these values 
are used for the MBM'B'j vertices (41), while for the 
baryon anomalous magnetic moments we use the exper- 
imental values instead of Ki (42). For the unobserved 
S° magnetic moment, we use the SU(3) flavor relation 
/i^o = (/is+ -l- fi-s-)/2, which is consistent with quark 
models. The transition magnetic coupling S*^ — >■ A7 
does not contribute in the isospin symmetric limit, be- 
cause this interaction changes the isospin of the excited 
baryon to 1. Now the total BB'-f coupling Vg, and the 
M BM'B'^ vertex F^- are given by the sum of the normal 
and anomalous contributions: 



2M, 



2M,, 



2Afn 



"^^'^'KQT.+gxJ-A.,^ , (45) 



4/2 



2M„ 



where /ie- is the observed magnetic moment of the 
baryon. 

Here we note that the magnetic interactions of the ex- 
cited baryons are obtained only from the baryonic dy- 
namics, since in our approach the spinless meson is bound 



the electromagnetic field tensor F^,^^ the charge matrix 
Q, the spin matrix 5*^, the SU(3) matrix of the baryon 
octet field B, and the chiral field = U = exp(i-\/2$//) 
where $ is the SU(3) matrix of the Nambu-Goldstone 
boson field. This interaction Lagrangian gives us spa- 
tial components of both the BB'j and the MBM'B'^ 



Here wc note that in Ref. [18] the additional contribution (Afeg = 
1) was introduced in order to take into account the normal mag- 
netic part of the gauged Weinberg- Tomozawa interaction. This is 
unnecessary in our calculation, where 6g = 0.82, because Eq. (37) 
does not produce the magnetic part. 
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by the baryon in s-wavc channel and the spatial compo- 
nent of the MM' J couplings do not contribute to the 
magnetic interactions in the Breit frame of the excited 
baryons. 

In order to study the internal structure of the A(1405) 
theoretically, we also consider the form factors probed 
with the baryonic and strangeness currents. For the bary- 
onic and strangeness current interactions, we replace the 
meson and baryon electric charges with the correspond- 
ing quantum numbers, namely, the baryon number and 
strangeness of the mesons and baryons: 



iM 



0, 



1, 



for the baryonic current and 



QiT = Qij = Qn = 0, 
Qk ^ Qk = Qt. = —1, 
Qk = 1, <3s = ^2, 



(46) 



(47) 



for the strangeness current. Then we consider the 
time components as the form factors for the baryonic 
and strangeness current interactions, as described below 
Eq. (6). Since the baryonic and strangeness form factors 
are calculated by the time component of the current, we 
do not need the counterparts of the matrices Xij and Yij 
given in Eq. (43), which contribute to the spatial compo- 
nent. 



E. Calculation of photon-coupled meson-baryon 
amplitudes 

Now we discuss the details of the calculation of the 
scattering amplitude of the MBj* — > M'B' process in 
the chiral unitary approach, in which the amplitude for 
the AIB — M'B' is given by multiple scattering of the 
meson and baryon. 

One of the most important issues wc take account of 
is the charge conservation in the calculation of the scat- 
tering amplitude for the AIB-f* — >■ M' B' process. This 



ensures the correct normalization of the form factor of 
the excited baryon, F^iQ^ = 0) = Qem, Fb{Q^ = 0) = 
B = 1, and Fs{Q'^ ~ 0) — S. Following the method pro- 
posed in Refs. [28, 48, 64], to calculate the form factors 
we take three relevant diagrams shown in Fig. 3, which 
contain the double-pole terms for the resonance states. 
Although charge conservation requires seven other di- 
agrams as shown in Fig. 4 for the general amplitude 
[64], these diagrams cannot contribute to the matrix 
elements at the resonance pole calculated by Eq. (19), 
since these terms have at most a single pole [28]. This 
means that, on the resonance pole, the charge conser- 
vation is maintained by only three diagrams shown in 
Fig. 3. Summing up the diagrams in Fig. 3, we obtain 
the relevant amplitude for the evaluation of the form fac- 
tors: 



rpfl _ rpfl 

7U — 7(l)ii 



(48) 

These contributions can be expressed by the combination 
of the BS amplitude and the elementary couplings dis- 
cussed before, according to the Feynman diagrams given 
in Fig. 3. In the Breit frame, in which the momenta of 
the photon and the A(1405) before photon coupling arc 

expressed as ~ (0, q) and P*" = ^y^s + q^/4, —q/2^, 

respectively, their explicit forms are written as: 



Kim = E T^k{V~s)D^^ {V~s; Q^)n, (v^), (49) 

k 

^7(2)..- = E T^k{Vs)D'^, (V^; Q^)nAV~s), (50) 

k 

(51) 



k,l 



where the vertex F^^ is given in Eq. (45) and the loop 
integrals with the photon couplings to the meson and 
baryon are given by. 



J {2tt [P - qiY ~ M^ + le [qi + qy - mf, + le^ '^"^ j _ 



k It - ml + ie 

2Mk 



D^^^ iVs; Q') = z / (2^)4 ,2 „ 4, + (P + g - q!r- Ml + ^e [^B. " ^ - ?1 + '?)] _ 



(52) 



(53) 



In the Breit frame s = P^P^ = (P + q)^{P + q)^ and related to the loop integral G given in (31) as, 
2P^(j',, = = > 0, and Dm, D-q, and F are functions 
of s/s and Q^. The function Dm (i's) at = is 

^^.(VJ;Q'-o) = gM.^, (54) 
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T- 



7(1) 




7(3) 



FIG. 3: Diagrams for the Ti^ which contain the double-pole terms of the excited baryon [28, 48, 64]. The shaded ellipses 
represent the BS amplitude. The dashed, solid and wiggly lines correspond to the ground state meson, the ground state baryon 
and the probe current, respectively. 



m 



T 



7(4) 



T 



7(5) 



7(6) 




T' 



7(7) 



7(8) 



7(9) 



T 



^7(10) 



FIG. 4: Supplemental diagrams for the charge conservation of the Ti^ in addition to the diagrams in Fig. 3 [64]. The shaded 
ellipses represent the BS amplitude. The dashed, solid and wiggly lines correspond to the ground state meson, the ground state 
baryon and the probe current, respectively. 



i^^,(Vi;Q' = o) = gB.^, (55) 

which can be easily proved by calculating the derivative 
of the loop integral G and using the vertices (35) and 
(44). Other analytic properties of these loop integrals, 
Dyi and Db, are discussed in Appendix B. ^ 

Now we introduce the finite size effects of the con- 
stituent hadrons in our scheme. The ground state mesons 
and baryons have spatial structures, while we have im- 
plicitly assumed that they are pointlike particles in the 
effective Lagrangian approach. Therefore we should in- 
clude the finite size effects of the constituent hadrons in a 
gauge invariant way. Here we simply multiply a common 
form factor (CFF) ^cff((3^) of the constituent hadrons 
to each photon vertex, 

DtiSV~s\ Q^) Dl^S^'-^ Q^)FcFFiQ^), (56) 
D^, (V^; Q') ^ Q')Fcff{Q'), (57) 

niiVs; Q2) ^ r^;(V^; Q^)FcMQ^)- (58) 



If one adopts the cut-off procedure for the meson-baryon loop 
integral G, one must calculate the photon-coupled meson-baryon 
loop integrals Dm and Db in a consistent way with the same 
cut-off. Form factors of A(1405) in such a cut-off procedure were 
evaluated in Ref. [30]. The two approaches, the dimensional 
regularization in our procudure and the cut-off procedure, would 
not make much differences in form factors except for the high 
momentum region compared to the cut-off scale. 



Note that the Fcff only depends on so it can be 
factorized out from the loop integrals of Eqs. (52) and 
(53). In this study we employ the dipole type form factor 
as, 

^CPP(Q^)^(^)'. (59) 

We take = 0.71 GcV^, which reproduces nucleon form 
factors well. This CFF corresponds to the hadron density 
^ exp(— Ar) with radial coordinate r of each hadron and 
mean squared radius (r^) = 12/A^ ~ 0.66 fm^. 

Here we comment on the normalization of the effective 
form factor F'^^ obtained on the real axis. As discussed 
in Appendix A, the effective form factor is correctly nor- 
malized so as to give the corresponding charge of the 
system at = independently of the energy yi, if 
we consider all of the diagrams shown in Figs. 3 and 4. 
Among them, the three double-pole diagrams in Fig. 3 
give dominant contributions at the energies close to the 
resonance pole position. Therefore, even if we take into 
account only the three diagrams in Eq. (48), deviation 
from the correct normalization for the effective form fac- 
tor is considered to be small. We also note that a part 
of the less singular terms in T^^j (hence F°^) is included 
in these three diagrams because off the pole position the 
nonresonant background in [see Eq. (11)] generates 
less singular parts in Eqs. (49)-(51). Here we choose the 
KN{I ~ 0)7* KN{I = 0) channel for the evaluation 
of F'^^ for the A(1405), which reduces the contribution 
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from the less singular terms thanks to the large coupling 
strength g^jvi demonstrated in Eq. (23). 

At last, let us show a relation among electric (Fe), 
baryonic (-Fb), and strangeness (Fg) form factors for the 
A(1405) based on isospin symmetry, in which we have a 
generalized Gell-Mann-Nishijima relation for the probe 
current. 

J^ = Jl+\[J^ + J^l (60) 

with the current for the third component of the isospin, 
. Taking the matrix clement for the A(1405) with 
7 = 0, we obtain'', 

F^m = \[F^{Q^)+Fs{Q^)l (61) 

where the matrix element of Jj vanishes for the 7 = 
state. From the relation (61), we have the relation for 
the spatial densities, 

PE(r) = i[pB(r) + ps(0], (62) 
with radial coordinate r as well. 

III. NUMERICAL RESULTS 

In this section, we discuss the internal structure of the 
resonant A(1405) state. We will show the numerical re- 
sults of the A(1405) form factors measured by the elec- 
tromagnetic, baryon number, and strangeness currents in 
momentum space. We will also show the spatial density 
distributions in coordinate space which are obtained by 
performing Fourier transformation of the form factors. 

The calculation of the form factor is performed in two 
ways; one is on the A(1405) pole position (Sec. Ill A and 
the other is to evaluate the effective form factor (21) on 
the real energy axis around the resonance energy region 
yfs - 1420 MeV (Sec. IIIB. On the A(1405) pole posi- 
tion, the internal structure of the resonance can be ob- 
tained in a way to keep charge conservation without non- 
resonant background contributions. The results, how- 
ever, may not be directly compared with experimental 
observables. The effective form factor on the real en- 
ergy axis, on the other hand, may be determined in ex- 
periments, but the obtained form factors have both the 
resonant and nonresonant contributions. 

For a reference of the electric size of the typical neu- 
tral baryon, we will compare the electric A(1405) form 
factor with an experimental fit of neutron electric form 
factor [65], 



At = 0, we obtain the usual Gell-Mann-Nishijima relation, 
Qem = {B + S)/2, for h=0 hadrons 



with a = 1.25, b = 18.3, A^ = 0.71 GeV^ the neu- 
tron mass A/n, and the neutron magnetic moment /in = 
— 1.913 /iN, where fi^ is the nuclear magneton. 

We also compare the results of magnetic, baryonic and 
strangeness form factors with a dipole type form factor, 

Fdipolc(g') = c X (^^2^) , (64) 

with A^ = 0.71 GeV^, which reproduces well the observed 
nucleon magnetic form factor [66]. The overall factor c 
will be adjusted to the normalization of the form factor 
of the A(1405) in question. 

A. Form factors on the resonance pole 

Here we discuss the internal structure of the A(1405) 
using the form factors obtained at the pole position in the 
complex energy plane. We evaluate the form factors of 
the higher A(1405) state, Z2, out of two A(1405) states, 
since this state gives the dominant contribution to the 
spectrum and is considered to be originated from the 
KN bound state. The form factors of the A(1405) at 
the resonance position arc obtained by Eq. (19) together 
with the amplitudes calculated in Eqs. (33) and (48). 

1. Electromagnetic, baryonic, and strangeness structures 

First of all, we show our results of the electric and 
magnetic form factors in Fig. 5 together with the em- 
pirical form factors of the neutron given in Eqs. (63) 
and (64). The normalization parameter is given by 
the real part of the magnetic moment of the A(1405), 
c = RelFuiQ^ = 0)]. Here, in order to see the finite 
size effects of the constituent hadrons, we also show re- 
sults without the common form factor (CFF) introduced 
in Eq. (59). The finite size effects make the magnitude 
of the form factors reduced, especially in the large 
region. Hereafter, we show only the results with CFF 
unless explicit mentionings. 

Now let us discuss the electromagnetic form factors of 
the resonant A(1405) shown in Fig. 5. The form factors 
-Fe and Fm contain the imaginary parts, since they are 
evaluated on the resonance pole position in the complex 
energy plane. However, we obtain the imaginary parts in 
smaller magnitude than the real parts. This is the conse- 
quence of the relatively small imaginary part of the pole 
position of Z2, since the form factors are real numbers 
in the limit of zero imaginary part of the pole position. 
For the charge neutral A(1405), deviation from zero in 
the electric form factor indicates that the A(1405) has a 
nontrivial charge distribution as seen in the neutron form 
factor. Comparing the real part of the A(1405) form fac- 
tor and the empirical neutron form factor, we find that 
the A(1405) form factor has larger magnitude than that 
of the neutron, especially at low Q^. This indicates that 
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FIG. 5: Electromagnetic form factors of the A(1405) state on the higher pole position Zi, together with the empirical form 
factors of the neutron. Left (right) panel shows the electric (magnetic) form factor _Fe (^ivi). The label "w/o OFF" represents 
the result without inclusion of the common form factor in Eq. (59) . The parameter c in the dipole form factor is chosen to be 
c — Re[_FM((5^ = 0)], the real part of the magnetic moment of the A(1405). 
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FIG. 6: Baryonic (Fb) and strangeness (Fs) form factors 
of the A(1405) state on the higher pole position Zi. The 
strangeness form factor is presented with the opposite sign 
for comparison. The parameter c in the dipole form factor is 
chosen to be c = 1. 

the spatial structure of the A(1405) is larger than the 
neutron. For the magnetic form factor of the A(1405), to 
which only the baryon components contribute as seen in 
Sec. II D, the real part of the result shows faster decreas- 
ing than the dipole fit to the nucleon. These results of 
i<E and i^M suggest the peculiar electromagnetic struc- 
ture of the A(1405) compared with the typical neutral 
baryon such as the neutron. 

Next we discuss the baryonic and strangeness form fac- 
tors of the A(1405) using the baryonic and strangeness 
currents as external probes. The calculation is done in 
the same way as the electric form factor using the bary- 
onic charge (46) and the strangeness charge (47) instead 
of the electric charge. The baryonic and strangeness form 



factors, F-Q and _Fs: are plotted in Fig. 6. For comparison, 
the strangeness form factor is presented with the opposite 
sign. Because of the baryon number and strangeness con- 
servation in our formulation the form factors are correctly 
normalized as ^b((3^ = 0) = 1 and i^s(Q2 = g) = -1, 
respectively. From Fig. 6 we find that the imaginary 
parts of the total baryonic and strangeness form factors 
are small compared with their real parts, as in the case 
of the electromagnetic form factors. We can also see 
that both the baryonic and strangeness form factors give 
steeper derivative at = compared with the nucleon 
form factor, and that the strangeness form factor shows 
faster decreasing than the baryonic one. Hence, the bary- 
onic and strangeness components also imply the peculiar 
structure of the resonant A(1405) among the ordinary 
low-lying hadrons. 

Now it is interesting to visualize the spatial structure 
of the resonant A(1405) in coordinate space as the den- 
sity distributions and the mean squared radii obtained 
from the form factors using Eqs. (8) and (9). We intro- 
duce a normalized density distribution P(r) = ^-kt^ p(r) 
with the density distribution p(r). With this definition, 
integrating P(r) from to oo, one gets the total charge 
(or magnetic moment) of the system: 

/•oo 

/ drV{r)=Q. (65) 
Jo 

In Fig. 7 we plot the electromagnetic density distribu- 
tions of the A(1405). As we can see from the left panel, 
the real part of the charge distribution Pe has large mag- 
nitude compared with that of the neutron. This may in- 
dicate smaller overlap between positively and negatively 
charged components in the A(1405) than in the neutron. 
It is also found that the charge distribution has positive 
values in the inner part (r < 1 fm) and negative val- 
ues in the outer part (r > 1 fm). Since the A(1405) 
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FIG. 7: Normalized distributions 47rr^p(r) of charge (Pe, left) and magnetic moment (Pm, right) densities of the A(1405) state 
on the higher pole position Z2 ■ Empirical charge distribution in the neutron is evaluated by Eq. (63) . Line denoted as "Dipole" 
in magnetic moment density is evaluated by Eq. (64) with c = Re[-FM(Q^ = 0)]. 
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FIG. 8: Baryonic (Pb) and strangeness (Ps) density distribu- 
tions of the A(1405) state on the higher pole position Z2. The 
strangeness density distribution is presented with the opposite 
sign for comparison. Typical density of nucleon is evaluated 
by Eq. (64) with c = 1. 



(Z2) strongly couples to the KN channel, the charge 
form factor of the A(1405) is expected to be dominated 
by the K~p component. Hence, our result implies that 
the lighter K~ surrounds the heavier p. From the mag- 
netic distribution Pm (Fig. 7, right) we can see the spa- 
tially larger structure of the A(1405) compared with the 
neutron. Since the magnetic density distribution is con- 
tributed mainly from the magnetic moment of the baryon 
in the A(1405), the magnetic structure reflects the bary- 
onic component. 

In the same way, the baryonic and strangeness den- 
sity distributions are plotted in Fig. 8 with the typical 
density distribution of the nucleon, which is evaluated 
by the Fourier transformation of the dipole form factor 



TABLE II: Electromagnetic (upper), baryonic and 
strangeness (lower) mean squared radii of the A(1405), 
{r'^)E, {r^)M, {r^)B, and {r^)s, on the higher resonance pole 
position Z2. 





-0.157 + 0.238i fm^ 




1.138 - 0.343i fm^ 




0.783 - 0.186j fm^ 




-1.097 + 0.6624 fm^ 



in Eq. (64) with c = 1. For comparison, the strangeness 
density distribution is shown with the opposite sign. The 
baryonic and strangeness density distributions clearly in- 
dicate the dominance of the real parts over the imaginary 
parts and spatially larger structure than the nucleon. 
In addition, the strangeness density distribution has a 
longer tail than the baryonic one. Since the A(1405) at 
pole position Z2 is dominated by the KN component, 
this means larger K distributions compared with N in- 
side the A(1405). 

Finally we evaluate the electromagnetic, baryonic, and 
strangeness mean squared radii, which are calculated us- 
ing Eq. (9) with the form factor. For the magnetic mean 
squared radius, we use the calculated magnetic moment 
Fm{Q^ = 0) = (0.17 — 0.05i)^N as the normalization of 
the mean squared radius. The results are shown in Ta- 
ble II. We find that the absolute value of the electric 
(magnetic) mean squared radius is |(r^)E| — 0.29 fm^ 
(|(''^)m| — 1-19 fm'^), which is about two times larger 
than that of the neutron ^ —0.12 fm'^ (~ 0.66 fm^). 
Also the absolute values of the baryonic and strangeness 
mean squared radii are larger than the typical size of nu- 
cleon. We also observe in Table II larger radius of the 
strangeness distribution than the baryonic one. This is 
due to the effect of the longer tail in strangeness density 
distribution compared with the baryonic one. Therefore, 
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these results support that the resonant A(1405) state has 
a spatially-extended structure compared with the typical 
baryon size < 1 fm. 

As a consequence, all of the results for the electromag- 
netic, baryonic, and strangeness structures show that the 
resonant A(1405) has a large size compared with typical 
hadrons. Furthermore, it is interesting to observe that 
the strangeness density distribution of the A(1405) has 
longer tail than the baryonic one. Since the A(1405) {Z2) 
is dominated by the KN component and K {N) carries 
the strangeness (baryon number), one can expect that 
such behaviors of the strangeness and baryonic distribu- 
tions are understood by the widely spread K distribution 
around N inside the A(1405). This expectation is sup- 
ported by the charge distribution, since it has positive 
values in the inner part whereas negative values in the 
outer part, which will be caused by the K^p component 
inside the A(1405). In the next subsection, we will clarify 
these detailed structure of the A(1405) by decomposing 
the form factors and density distributions into the con- 
tribution from each meson-baryon component. 



2. Contribution from each meson-baryon component 

In order to discuss the inner structure of the A(1405) 
resonance from the theoretical point of view, it is inter- 
esting to decompose the form factors into the contribu- 
tion from each meson-baryon state to which the external 
current couples. This decomposition can be done by cal- 
culating the amplitude T^^^^ -hT^^j) (^9) ^^'^ (^0) 
without the summation of the intermediate channel k, so 
that the total form factor is obtained by summing up 
all the components. In the decomposition we do not in- 
clude the contribution from T^^^y which is not important 
for the study of the spatial size of the resonance, since 
this is a contact interaction and gives an almost trivial 
momentum dependence from the form factor of the con- 
stituent hadrons^. Here we mainly discuss the KN and 
ttS components in the form factors, since the KN and ttS 
intermediate states are the dominant contributions in the 
A(1405) resonance as we have seen in Table I, and we will 
check that the -qK and KE. channels and the contact term 
coming from T!.y(3) give negligibly small contributions to 
the form factors. 

We first discuss the contributions to the electromag- 
netic form factors from the KN and ttE intermediate 
states. In the upper panel of Fig. 9, we show the electric 
form factor coming from the K~p and 7r+E~ states. As 
one can see from the figure, the electric form factor from 
the K~p state reproduces almost the total form factor 
given in Fig. 5. The reason is as follows; the neutral 
hadrons do not contribute to the electric form factor. 
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For the electric current, ^^(3) gives no contributions in case of 
neutral resonances. 
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FIG. 9: Meson-baryon components of the electromagnetic 
form factors of the A(1405) state on the higher pole posi- 
tion Z2- Components of the elecric (magnetic) form factor 
are shown in the upper panel (middle and lower panels). 



and the A(1405) resonance (Z2) has a tiny coupling to 
the i^+S" channel. Then, the sum of the contributions 
from K~p and tt^S^ components dominate the electric 
form factor. However, since we are working in the isospin 
symmetric limit, the tt^S^ states give exactly same con- 
tributions in magnitude for the isospin 7 = resonance 
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TABLE III: Values of the baryonic and strangeness form fac- 
tors at = 0. The separated contributions from each meson- 
baryon channel and contact term are also listed. Each meson- 
baryon channel is given in the isospin basis; for instance, KN 
represents the sum of the contributions from K~ p and K'^n. 



Component 


Fb(0) = -i^s(O) 


total 


1 


KN 


0.994 + 0.048i 


ttE 


-0.047 - 0.151i 


77A 


0.052 + 0.012i 


K~. 


-0.002 + 0.002i 


Contact 


0.002 + 0.089i 



with the opposite sign due to their electric charges. Thus, 
the sum of the tt^E^ states does not contribute to the 
total electric form factor for the A(1405) resonance and 
the K~j) state reproduces almost the total electric form 
factor. 

In the -K^Yr contribution, we observe the sudden in- 
crease of the real part and the peak structure of the imag- 
inary part seen at ~ 0.2 GeV^, which comes from the 
analytic properties of the loop function in which the 
external current attaches to the pion propagator. We will 
discuss the details in next subsection and Appendix B. 

The separated contributions to the magnetic form fac- 
tor are shown in the middle and lower panels of Fig. 9. 
For the magnetic contribution, the meson pole term T^^j^) 
shown in Fig. 3 does not contribute and we are left with 
the baryon pole term ^"^(2) {^"^^ t^i^ contact term ^'^(3) , 
which contribution is not included in Fig. 9. Here we 
show the form factors with the K^n^ and 7r+S~ 

intermediate channels. We find that the KN contribu- 
tion is substantially larger than ttE. Here we also find 
that the K^j) contribution has the opposite sign to the 
K^n and they largely cancel each other. This is because 
only the isosinglet component can couple to the A(1405) 
due to the isospin symmetry and the isosinglet magnetic 
moment of nucleon is known to be very small. This result 
indicates that the magnetic form factor of the A(1405) is 
composed mainly by nucleons in KN dynamics in which, 
however, large cancellation between K"^ and K^n com- 
ponents takes place. 

Next we discuss the contributions to the baryonic and 
strangeness form factors from each meson-baryon chan- 
nel. Here, instead of plotting the form factors with re- 
spect to Q^, we show in Table HI the values of the bary- 
onic and strangeness form factors at = with indi- 
vidual contributions from the meson-baryon channels to 
the A(1405), which corresponds to the each channel con- 
tribution to the baryon number and strangeness of the 
system. Here we have the Gell-Mann-Nishijima relation 
Fb = — Fg at = for each meson-baryon channel. 
The decomposition to each meson-baryon channel implies 
that the KN{I = 0) channel gives more than 90% of the 
total baryonic and strangeness charges, whereas the ttE, 
TjA, /"TS, and contact-term components in / = chan- 



nel are negligibly small. The magnitude of the charge of 
each component is determined by the coupling strength 
given in Table I and the derivative of the loop function 
which contributes to the form factor through Eqs. (54) 
and (55). Thus, due to the large coupling strength .g^fjv 
the KN channel dominates the structure of the A(1405). 

Now we show the meson-baryon components of the 
electromagnetic density distributions, Pe and Pm, which 
are obtained by the Fourier transformation of the cor- 
responding meson-baryon components of the form fac- 
tors. The results are shown in Fig. 10. Note that Pe 
is plotted up to r = 10 fm instead of 5 fm. We find 
again the dominance of the KN component in both Pe 
and Pm- Hence, the negative (positive) charge distri- 
bution of the A(1405) in Fig. 7 is now understood as 
the lighter K~ (heavier p) existence in outside (inside) 
region. It is a more interesting finding that the tt+S^ 
(equivalently 7r~E+ with the opposite sign) component 
of the electric density distribution shows a characteris- 
tic behavior of dumping oscillation. As we will discuss in 
next subsection and Appendix B, the oscillating behavior 
can be interpreted as the decay of the system into the ttS 
channels through the photon coupling to the intermedi- 
ate meson. Although this oscillation is interesting from 
the theoretical point of view, it does not contribute to 
the total density distribution due to the cancellation of 
7r+E~ and 7r~E"'' components. For the magnetic density 
distribution, we again observe large cancellation between 
K~p and K^n components in the isosinglet A(1405). 

Then, we decompose the baryonic and strangeness den- 
sity distributions into the different meson-baryon contri- 
butions in Fig. 11. In this figure, KN represents the sum 
of the contributions from K~p and K^n, whereas ttE rep- 
resents the sum of the contributions from tt^E", tt+E", 
and 7r~S+. For the ttE channel, we plot the baryonic 
density which is equivalent to the strangeness ttE density 
with the opposite sign. This is clear from the charges in 
Eqs. (46) and (47) where both baryonic and strangeness 
currents probe the E component with the same strength 
and the opposite sign in the ttE intermediate state. As 
one can see from the Fig. 11, the KN component where 
the baryonic (strangeness) current probes the N (K) dis- 
tribution dominates both the baryonic and strangeness 
density distribution in Fig. 8. Therefore, the longer tail of 
the strangeness distribution in Fig. 8 is now understood 
as the larger distribution of the K component than that 
of the baryon number distribution generated by the N. 
This is consistent with the electric density distribution, 
that the lighter K^ locates outside the p, which should 
also be the case for the K'^ and n through the isospin 
symmetry. 

From the decomposition of the form factors and density 
distributions into meson-baryon components, it is con- 
firmed that the resonant A(1405) is indeed dominated 
by the KN{I = 0) component, giving more than 90% 
of the total baryonic and strangeness charges, and it is 
found that the A(1405) is dominantly composed of the 
K in the outside region around the nucleon, with a large 



16 



0.15 



0.1 



S 005 



-0.05 



-0.1 



Re. if p 
Im.K p 
Re. 71*2: 
Im. 71*2: 



\ 'I 



\: / 



4 6 

r [fm] 



10 



Re. if p 

Im.ifn 

Re.X°n 

Im.X°7i 




FIG. 10: Meson-baryon components of the distributions 
47rr^p(r) of charge (Pe, upper) and magnetic moment (Pm, 
middle and lower) densities of the A(1405) state on the higher 
pole position Z^- 



size compared with typical hadrons. We also find that 
the magnetic moment of the A(1405) is composed mainly 
by nucleons in KN dynamics in which large cancellation 
between K~'p and K^n components takes place. In ad- 
dition, we observe that the tt+E" component has the 
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FIG. 11: Meson-baryon components of baryonic (Pb) and 
strangeness (Pg) density distributions of the A(1405) state 
on the higher pole position Z2- The strangeness density dis- 
tribution is presented with opposite sign for comparison. Here 
KN represents the sum of the contributions from K~ 'p and 
K^n^ whereas ttE represents the sum of the contributions from 
TT^E", TT+E" and 7r-E+. 



escaping oscillation in the charge distribution (see also 
next subsection), which is although not observed in the 
total charge due to the tt^S^ cancellation. In our results 
all of the results for the electromagnetic, baryonic, and 
strangeness structure of the A(1405) are consistent with 
each other. 



3. Escaping oscillation in decay channel 

In the previous subsection we found an interesting be- 
havior of the 7r+S~ escaping oscillation in the electric 
density distribution. Here we discuss the escaping os- 
cillation found in the ttE channel, or decay channel in 
general. 

From a viewpoint of the Fourier transformation, what 
makes oscillation behavior in the tt~^Y,~ electric density 
distribution is the peak structure at ~ 0.2 GeV^ in 
the electric form factor in the tt+S^ channel. Namely, 
the Fourier transformation of the form factor picks up 
much contributions from the peak structure, making 
large magnitude in the corresponding wave number for 
the density distribution, which is seen as the oscillation 
component. 

Then, let us discuss the origin of the peak structure in 
the form factor. An important point is that the photon 
coupled loop integral of the intermediate channel k, Dm^ 
(Z?Bfc ) , contains a divergent point at 



Mi 



(66) 



with real energy ^/s, defined in Eq. (29), and the 
baryon (meson) mass {m-k) hi the intermediate 
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state. This divergence point corresponds to the t- 
channel threshold. The detailed discussion is given in 
Appendix B. It should be emphasized that this singular- 
ity (66) can be reached only in case that the energy is 
above the threshold of the k channel, ^/s > Mk + rtik , so 
that iq^s > 0. Such a singularity at certain gener- 
ates peak structure in the form factor with complex en- 
ergy through the analytic continuation ^/s — > z. Hence, 
that the resonance energy is above the threshold is es- 
sential to the peak structure in the form factors, and 
the oscillation behavior of the density distributions can 
be interpeted as the decay of the system into the open 
channels through the photon coupling to the intermedi- 
ate state, with kicked meson and baryon to the on-shell 
by the photon coupling. In the present case, since only 
the ttS channel is open for the A(1405) decay, the peak 
structure appears only in the ttE channel. 

From the above discussion, it is obvious that the peak 
structures in the form factors and the oscillation behav- 
iors in the density distributions should appear in decay 
channels for resonance states in the meson-baryon pic- 
ture. Such structures are, however, eventually not ob- 
served in the total electric form factor and density dis- 
tribution for the resonant A(1405). This is because, as 
we have mentioned, the tt^E^ components largely can- 
cel each other for the electric structure due to the isospin 
symmetry. Hence, if we would observe a excited state for 
which its decaying channel contributes to the form fac- 
tors without cancellation, we could observe the escaping 
oscillation in the total density distribution. 

At last, we comment on that the magnetic, baryonic, 
and strangeness density distributions do not show the 
(visible) oscillation behavior even in the decay chan- 
nel, ttE. This is because the magnetic, baryonic, and 
strangeness currents couple to E rather than pion in the 
ttE channel; the current coupling to the E propagator 
also provides the oscillation behavior in the electric as 
well as the magnetic, baryonic, and strangeness density 
distributions. However, Eq. (66) indicates that, due to 
the small pion mass to^, the E-current coupling makes 
the peak structure in the form factor with very high 
value [Q2 _ 10 GeV^ for the ttE channel in the A(1405) 
{Z2)], thus such a high coupling should be strongly 
suppressed by form factors of the constituent hadrons. 
As a consequence, the oscillation contributions from the 
E-current coupling are numerically small and not visible 
in Figs. 10 and 11. 



B. Effective form factors on the real energy axis 

Now we evaluate the effective form factors of the 
A(1405) on the real energy axis defined in Eq. (21), in 
order to study how the form factors obtained at the res- 
onance position are seen on the real energy axis. These 
effective form factors on the real axis will provide the 
quantities which can be compared with the experimen- 
tal observations, in contrast to the form factors obtained 



Re. part 

Im. part 





1 






CN 


0.8 









0.6 


II 






0.4 








0.2 









-0.2 
1 



390 1400 1410 1420 1430 1440 



s"^ [MeV] 

FIG. 12: Effective baryon number of the resonance, F^i^Q^ = 
0), as a function of meson-baryon center-of-mass energy y/s. 
Calculations are performed in the KN{I = 0)7* — >■ KN{I = 
0) process on the real energy axis. 



in the complex energy plane. Here we comment that the 
mean squared radii are not well defined for resonances on 
the real energy axis due to insufficient fall-off of the den- 
sities at large r, in contract to the case on the resonance 
pole position. 

As discussed in Sec. HE, the form factors calculated 
with three diagrams shown in Fig. 3 keep correct nor- 
malizations only on the resonance pole position. Off the 
resonance pole position the additional diagrams shown in 
Fig. 4 are necessary to keep the correct normalizations. 
Nevertheless, since the double-pole diagrams dominate 
the photon-coupled amplitude at the resonance energies 
if the pole position is not far from the real axis, the con- 
tributions from the supplemental diagrams in Fig. 4 to 
the form factors may be negligible. 

This can be checked by the calculation of the effec- 
tive baryon number F^^{Q^ = 0) for the A(1405) state 
with the diagrams in Fig. 3 on the real axis. Here we 
choose the KN — )- KN channel in / = with the bary- 
onic current, so that the less singular contributions are 
suppressed through the large coupling constant g^j^ [see 
Eq. (23)]. In Fig. 12, we plot Fg^{Q^ = 0) for the A(1405) 
on the real energy axis as a function of -y/s. which is eval- 
uated with the three diagrams shown in Fig. 3. If the 
amplitude maintains the conservation law by taking ac- 
count of all the diagrams, F^^{Q^ = 0) should be pure 
real and unity independently of ^/s. Figure 12 indicates 
that the real part of F^{Q^ = 0) calculated with the 
three diagrams is close to one in the energy A(1405). 
This means that the three diagrams dominate the ampli- 
tude of the kN[I = 0)7* KN{I = 0) process, and 
that contributions coming from the neglected diagrams 
are less than 10%. To study the energy dependence, we 
evaluate the form factors at three meson-baryon center- 
of-mass energies ^/s = 1410, 1420, and 1430 MeV. 

The results of the real part of the effective electromag- 
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FIG. 13: Real part of the effective electromagnetic form factors (F^ and ) on the real energy axis, together with the 
empirical form factors of the neutron. Calculations are performed with the center-of-mass energy -^s = 1410, 1420, and 
1430 MeV. The parameter c in the dipole form factor is chosen to be c = Ke[F^{Q'^ — 0; ^/s = 1430 MeV)], the real part of 
the magnetic moment at ^/s — 1430 MeV. 



netic form factors on the real energy axis are shown in 
Fig. 13. We find that the efi^ective electromagnetic form 
factors for ^/s = 1420 MeV are qualitatively very similar 
with the form factors evaluated at the pole position in 
Fig. 5. We also find in Fig. 13 that the effective electro- 
magnetic form factors have mild energy dependence even 
in energies close to the resonance position. This is be- 
cause, off the resonance pole position, less singular terms 
also contribute to the effective form factors and brings 
energy dependence to them. We also see that, in the en- 
ergies closer to the KN threshold, both the F^^ and 
have steeper change at low Q^. 

Now let us discuss the electromagnetic density distri- 
butions Pf and Pm obtained from Ff and F^. Wc plot 
the real parts of Pf and Pm in Fig. 14 with = 1410, 
1420, and 1430 MeV, in comparison with the empirical 
neutron density distribution. The important point is that 
the effective distribution in the real axis also catches the 
properties found in the density distribution on the pole 
position, that is, spatially larger structure than the neu- 
tron with the outward negative charge of K~ and the 
inward positive charge of p. As the energy -^s increases, 
the observed distribution becomes effectively wider. 

Finally we summarize the numerical analyses of the 
internal structure of the resonant A(1405) state, which 
is done by evaluating two observables; one is the reso- 
nance form factor (19) on the resonance pole position 
(see Sec. Ill A, and the other is the effective form fac- 
tor (21) on the real energy axis (Sec. IIIB. Although 
there is a quantitative difference between these two anal- 
yses such as the energy ^/s dependence of the effective 
form factors, it is found that the several peculiar fea- 
tures of the resonance form factor in the complex plane 
are mostly maintained in the effective form factor on the 
real axis. This means that the properties of the A(1405) 
defined in the complex plane may be within our reach of 



the experimental searches achieved by the real energies. 
This consequence comes from the small imaginary part 
of the pole Z2 of the A(1405). 

From these analyses we have drawn the following con- 
clusions. The A(1405) has softer form factors than those 
of neutron. Consequently the mean squared radii are 
larger than the neutron radii which can be regarded as 
a typical baryon, as found in Rcf. [28]. Through the de- 
composition into meson-baryon channels and the analysis 
of the form factors with different probe currents, the in- 
ternal structure of the resonant A(1405) is found to be 
dominated by the KN component, with the nucleon in 
the center being surrounded by the antikaon. 



IV. DISCUSSIONS 

We have considered the form factors and density dis- 
tributions of the resonant A(1405) state. Because the 
A(1405) resonance has the finite decay width, the ob- 
tained form factors are complex numbers. Although we 
have deduced the structure of the A(1405) mainly from 
the behavior of the real part of the form factors, the in- 
terpretation may not be as straightforward as the stable 
particle. 

The A(1405) has been considered as a quasi-bound 
state of the KN system^ having the ttS decay chan- 
nel [1, 2]. This picture is also supported by the chiral 
unitary approach, where the KN bound state is gener- 
ated only by the attractive KN interaction and chan- 
nel coupling to ttS provides the A(1405) with the decay 



In this section KN represents KN{I = 0) channel in the isospin 
basis. 
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FIG. 14: Real part of the effective charge (Pg , left) and magnetic moment (P^ , right) density distributions on the real 
energy axis. Calculations are performed with the center-of-mass energy -^s = 1410, 1420, and 1430 MeV. Empirical charge 
distribution in the neutron is evaluated by Eq. (63). Line denoted as "Dipole" in magnetic moment density is evaluated by 
Eq. (64) with c = Re[F^{Q'^ = GeV^; ^ = 1430 MeV)]. 



width [25]. In the present study, we confirmed this pic- 
ture by the decomposition of the baryonic charge in the 
previous section. Therefore, it is instructive to under- 
stand the structure of the KN bound state, by switching 
off the couplings of KN to other channels. 

In the case of the bound sate, the form factors are 
obtained as the real numbers for all , hence we can in- 
terpret the form factors as physical quantities. Further- 
more, the single-channel model allows us to investigate 
the internal structure of a dynamically generated bound 
state. For instance, we expect that the spatial size will be 
larger for the system with smaller binding energy. This 
does not trivially follow from the calculation of the form 
factor and is worth examining in the present framework. 
Thus, we examine the structure of the KN bound state 
by changing the model parameters. 

In Sec. IV A, we will study the structure of the 
KN bound state with the natural subtraction constant 
a^fjv — —1-95, which is obtained so as to exclude ex- 
plicit pole contributions from the loop function [9] . This 
ensures that the bound state has the pure molecule struc- 
ture. Then we will study the structure of the dynamically 
generated KN bound states with different binding ener- 
gies in Sec. IV B. The interaction strength is adjusted to 
fix the binding energy within the natural condition. We 
also discuss the meson-nuclcon bound state with differ- 
ent meson masses instead of the physical kaon mass, in 
order to see the effect of meson masses to the bound state 
in Sec. IV C, where the spatial size of the bound state is 
kept fixed. 



A. Structure of the KN bound state 

In this subsection, we consider the ^A'' single chan- 
nel model with a bound state. We use the Weinberg- 



Tomozawa interaction (25) for the KN channel and the 
subtraction constant = —1.95 with the regulariza- 
tion scale /iicg = 630 MeV, which is obtained by the con- 
dition G f^pf{AIisi) = [9]. With this subtraction con- 
stant, solving the ^A^ single-channel scattering equa- 
tion (33) we obtain the KN bound state at 1424 MeV 
with the coupling strength to the KN [sec Eq. (10)] as 
gjifq = 2.17. This bound state is generated only by the 
attractive A'A^ interaction and has a binding energy of 
10 MeV. 

Here we first study the structures of the K and A^ 
components in the A'A^ bound state. Since the K {N) 
has strangeness —1 (0) and baryon number (1), we can 
observe K and A^ distributions in the bound state by 
using the strangeness probe current with the opposite 
sign and the baryonic current, respectively, which leads 
to the relations for the form factors as, 

Fk{Q^)^-Fs{Q^), Fn{Q')^Fb{Q'), (67) 

where -f/f (jv) represents the K (N) form factors probing 
the K (N) component in the bound state. We also study 
the electric structure of the bound state, which is related 
to the K and A^ form factors via the generalized Gell- 
Mann-Nishijima relation (61). Here we do not consider 
magnetic component of the bound state. In Fig. 15, the 
K, N, and electric form factors are shown together with 
the empirical form factors of the neutron (63) and (64). 
Comparing with the previous results for the resonant 
A(1405), we find that all of the electric, A^, and K form 
factors of the KN bound state show similar characteris- 
tic behaviors with the real part of the electric, baryonic, 
and opposite-sign-strangeness form factors of the reso- 
nant A(1405), respectively. For the K and A^ form factors 
Fk,n^ we observe the steeper derivative at = com- 
pared with the dipole form factor and the faster decreas- 
ing K form factor than the A^ one, which reflect a salient 
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FIG. 15: The K, N (left) and electric (right) form factors of the KN bound state with mass 1424 MeV, together with the 
empirical form factors of the neutron. The parameter c in the dipole form factor is chosen to be c = 1. 



TABLE IV: The K, N, and electric mean squared radii, 
(r^)/^, {r'^)N, and (r^)E, and mean squared distance between 
K and N, {x^)kn, of KN bound state with mass 1424 MeV. 





{r2)iv [fm^] 




{x^)iiM [fm^] 


1.878 


0.998 


-0.440 


2.848 



structure of the bound state, as the resonant A(1405) dis- 
cussed before. Also the electric form factor exhibits 
large enhancement at smaller Q region (< 0.1 GeV^) 
and slowly decrease above > 0.2 GeV , as the real 
part of the form factor of the resonant A (1405). 

Using the form factors we calculate the density distri- 
butions through the Fourier transformation. In Fig. 16, 
the K, N, and charge density distributions from the 
form factors arc shown with the normalization P(r) = 
Airr'^ p{r). As we can sec from the figure, the behaviors of 
the density distributions arc similar to those for the res- 
onant A(1405) in Figs. 7 and 8, thanks to the large KN 
coupling gj^j^ of the A(1405); the charge distribution has 
positive values in the inner part whereas the negative val- 
ues in the outer part, and the K distribution has longer 
tail than the N one. 

The results of the K, N, and electric mean squared 
radii evaluated from the form factors are displayed in 
Table IV. The results of the K and N mean squared 
radii, (r^)^ = 1.878 fm^ and {r^)N = 0.998 fm^ re- 
spectively, indicate that both the K and N distributions 
spread compared with the typical nucleon size in the KN 
bound state, and the K has larger distribution than the 
N inside the bound state. The electric mean squared 
radius (r^)E is —0.440 fm^, whose absolute value is four 
times larger than that of the neutron ~ —0.12 fm^. In- 
deed, the charge density distribution Pe clearly shows the 
larger structure of the KN bound state than the typical 
neutral hadron. 

Next, it is also interesting to evaluate the mean squared 



distance between K and N from the form factors in our 
approach. For this purpose we take the nonrelativistic 
limit and treat the constituent hadrons as point particles, 
neglecting the common form factor effects [see Eq. (59)]. 
In the nonrelativistic limit, both the K and N form fac- 
tors Fff and Fjv are determined from the form factor 
for the relative motion in the two-body system i^rei with 
appropriate scale transformations as. 




M, 



N 



(68) 



(69) 



due to the kinematics of the system. Therefore, the 
mean squared distance between K and TV, {x'^)kn ~ 
6 dFrci/(i(3^|Q2=0i can be determined from both the K 
and N form factors Fj^{Q^) and F^iQ"^) with appropri- 
ate coefficients. The results are given by. 



N 



N 



2.849 fm^ 



X {r^)N = 2.846 fm^ 



(70) 



(71) 



for the mean squared distance between K and N. The 
two values in Eqs. (70) and (71) show very good agree- 
ment with each other. A small difference between two 
values in Eqs. (70) and (71) is expected to come from 
the field theoretical evaluation of the form factor and the 
relativistic correction with respect to the binding energy. 
In this study we define the mean squared distance as the 
average of the values evaluated from the K and radii. 
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FIG. 16: The K, N (left) and charge (right) density distributions of the KN bound state with mass 1424 MeV. Typical density 
of nucleon is evaluated by Eq. (64) with c = 1. Empirical charge distribution in the neutron is evaluated by Eq. (63). 



and we obtain the mean squared distance between K and 
N as 



(x2)^jv = 2.848 



(72) 



for the KN bound state with binding energy 10 MeV. 



B. KN system with different binding energies 

Here we discuss the KN bound state with different 
binding energies, in order to see the dependence of the 
structure on the binding energy. This can be achieved by 
replacing the interaction strength C = 3 for KN{I = 0) 
in Eq. (25) with Cb representing an interaction strength 
for KN to generate a bound state with a specified bind- 
ing energy Be, and we keep the subtraction constant 
a^N — —1-95 in order to exclude the explicit pole con- 
tribution. In Fig. 17 we plot interaction strength Cb as 
a function of the binding energy Be- 

We dynamically generate a KN bound system with 
different binding energies Be and show the electric, K, 
and N mean squared radii as functions of the binding 
energy in Fig. 18. It is obvious that the distribution of 
the constituent hadrons in the KN bound state spreads 
if the binding energy decreases, in accordance with our 
expectation from quantum mechanics. In addition, the 
mean squared radii are much sensitive to the binding en- 
ergy in the near-threshold region (B^ ^ 10 MeV), which 
indicates that for the shallow bound state the binding 
energy is a key quantity for the spatial structure of the 
bound state. Wc also find that the K distribution in the 
KN bound system is more sensitive to the binding en- 
ergy in this region, since the K is lighter than the N. 
In the large binding energy region, the K and N mean 
squared radii asymptotically goes to finite values. This 
behavior in the large binding energy region is caused by 
that the distributions of both K and N shrink to the fi- 
nite K and N radii. In contrast, the electric radius goes 
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FIG. 17: Interaction strength needed to generate a bound 
state with different binding energies, Cb- In the figure, 
Mbound represents the mass of the bound state. 



to zero, because K and N get close to each other, which 
leads to almost zero electric mean squared radius. 

Next let us make a simple comparison of our results 
of mean squared radii, which is based on the field the- 
ory, with that of the nonrelativistic wave function of the 
two-body bound state, in order to check the consistency 
of these methods. Outside the interaction range i?int the 
nonrelativistic wave function for the bound state takes 
an asymptotic form ^ exp {^ — ^2^Bex) jx in the rela- 
tive coordinate with the distance x and the reduced mass 
[I = MNmjf/{Mpf + nij^) = 324 MeV. Here we adopt a 
simple form having the correct asymptotic behavior, 

^(x) = (const.) X --H-V^^) ^ (73) 

X 

as a typical wave function for comparison. Using this 
wave function ip{x), one can evaluate the mean squared 
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FIG. 18; Mean squared radii of KN bound state with different 
binding energies. In the figure, Afbound represents the mass 
of the bound state. 



FIG. 20: Interaction strength Cb as a function of the meson 
mass m, in order to generate a bound state with the condi- 
tion (75). 
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tides, with that obtained from the nonrelativistic wave 
function (74), (a;^)NR. As we can sec from the figure, the 
difference between {x'^)xn ^^'^ (a^^)NR is quite small for 
every binding energy. This is an interesting result, be- 
cause it is nontrivial that the mean squared distance be- 
tween K and in our approach, where the bound state 
is dynamically generated and the size is probed by ex- 
ternal currents within the framework of the field theory, 
shows very similar value to that from the nonrelativistic 
wave function ip. 
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FIG. 19: Mean squared distance between K and A'' in the 
bound system in our approach (solid line), together with that 
obtained from the nonrelativistic wave function ip(x) (dashed- 
dotted line, denoted as "NR"). In the figure, Afbound repre- 
sents the mass of the bound state. 



distance of the two particles in the bound state as. 



dx Attx'^ \ iIj{x)\' 



1 



NR 



dxA-Kx^ \i^{x)\ 



2 4//Be 



(74) 



Although it is derived from the specific wave func- 
tion (73), this relation holds model- independently for the 
state with small binding energy, where the inner part 
{x < i?int) of the wave function is irrelevant and the 
asymptotic form of the wave function dominates the in- 
tegration in Eq. (74). Now in Fig. 19 we compare our 
result of the mean squared distance between K and N, 
(x^)^jYi which is evaluated by using method developed 
in the previous subsection for a bound state of point par- 



In this subsection, we vary the meson mass from the 
physical kaon mass, in order to see the effect of the me- 
son mass to the meson-nucleon bound state. Here we con- 
sider the meson (m~ , m'^) having same quantum numbers 
as the antikaon {K~, K'^) except for the mass. 

During the calculation, we fix the spatial structure of 
the bound system so that the meson-nucleon distance is 
to be in hadronic interaction range. For this purpose, we 
impose the condition. 



^Be = const. 



(75) 



where _Be is the binding energy of the system and /i is 
the reduced mass, /i = niMN/{m + Mjv)^*^. The in- 
teraction strenght Cb is adjusted by the binding energy 
through (75) and the natural subtraction constant de- 
termined by Gxjq{y/s = Mm] rrij^ — >■ m) = [9]. We 
start with the bound state with the physical kaon mass 



For the study of the different meson masses, we denote the meson 
mass by m. 
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FIG. 21: Charge density distribution of meson-nucleon bound 
state with different meson masses, 200, 500, 1000, 1500, and 
2000 MeV. 



TO = = 495.67 McV which leads Be = 10 MeV as 
studied in section IV A. In Fig. 20 we plot interaction 
strength Cb as a function of the meson mass. 

For the meson-nucleon bound state, the electric probe 
is useful to study the structure of the system, since the 
electric current probes the negatively charged m~ and 
positively charged p in the bound state, hence we can ob- 
serve relative distribution of the meson and nucleon in the 
system as well as the spatial size of the system. In Fig. 21, 
we show the charge density distribution of the meson- 
nucleon bound state for different meson masses. As we 
have discussed for the KN bound state, in the meson- 
nucleon bound state, concentration of positive (negative) 
charge density corresponds to the region dominated by 
the p (to~). For the small meson mass region (m = 200, 
500 MeV), the meson distributes at larger distance than 
the N, as the KN bound state in the preceding subsec- 
tions. If the meson mass is compatible with the nucleon 
mass (m = 1000 MeV) the charge distribution is almost 
flat around zero, which indicates that meson and TV have 
almost same distribution, and cancel each other in the 
charge distribution. For the meson mass larger than the 
nucleon (to = 1500, 2000 MeV), the meson and N distri- 
butions interchange their role, where the lighter nucleon 
goes outside and the heavier meson stays inside. Such 
a behavior is consistent with our expectation from the 
classical picture and indicates that we correctly observe 
internal structure of the bound state with respect to the 
mass ratio of the constituent hadrons. 

Here we note that the difference of the structures be- 
tween the bound states of KN and heavy meson-nucleon, 
if existed, might be realized as a difference of the be- 
haviors of the bound states in nuclear medium. For the 
KN (heavy meson-nucleon) bound state, N has smaller 
(larger) density distribution in the system, which im- 
plies that N has larger (smaller) momentum distribution 
due to the uncertainty principle. Therefore, in nuclear 



FIG. 22: Electric, baryonic and strangeness mean squared 
radii of meson-nucleon bound state with different meson 
masses. We also plot mean squared distance between meson 
and nucleon, (x^). 



medium, where the momentum of N is filled up to the 
Fermi momentum, the existence of heavy meson-nucleon 
bound state might be suppressed even in the lower den- 
sity compared with KN. 

It is also illustrative to plot the various mean squared 
radii and the mean squared distance between meson and 
nucleon as functions of the meson mass (Fig. 22. The 
mean squared distance is almost independent of the me- 
son mass, as a consequence of the fixed fiB-E- This sug- 
gests that the spatial size of the bound state is indeed 
constrained by this condition. The suppression of the 
mean squared distance in the smaller meson mass re- 
gion (to ~ 200 MeV) is caused by the condition that the 
binding energy, which is restricted by Eq. (75), is not neg- 
ligible compared with the meson mass. In addition, the 
electric, mesonic, and nucleon mean squared radii behave 
consistently with our interpretation discussed above. The 
sign change of the electric radius at m = Mjv reflects the 
inversion of the spatial distribution of the meson and N. 

At last we comment on the relation between the me- 
son mass m and the interaction strength Cb- The re- 
sult in Fig. 20 shows that in the small meson mass re- 
gion a strongly attractive interaction is required to keep 
the light meson binding with the nucleon, whereas the 
heavier meson can be bound as the molecular type (with 
(x^) ~ 3 fm^) by the moderate attraction, in accordance 
with the classical expectation. From the group theoret- 
ical point of view, the coupling strength of the scatter- 
ing of the octet baryon and the octet meson is at most 
C = 6 [67, 68]. The resuh in Fig. 20 suggests that if 
the meson were as light as 200 MeV, we would need the 
coupling strength C ^ 15 to generate a bound state with 
the size (x^) ~ 3 fm^. This means that the pion is too 
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light to generate a bound state of the molecular type"'^"'^. 
On the other hand, the sufficiently heavy meson mass 
(to > 400 MeV) is essential to generate a meson-nucleon 
bound state with spatial size of a hadronic molecular 
state by using the chiral interaction. 



V. CONCLUSION 

We have studied the internal structure of the reso- 
nant A(1405) state with electromagnetic, baryonic, and 
strangeness probes. The resonant A(1405) state has been 
described by the meson-baryon interaction picture based 
on chiral dynamics. The form factors and the density 
distributions of the A(1405) have been defined carefully 
as an extension of ordinary stable particles and been di- 
rectly evaluated from the meson-baryon scattering am- 
plitude, paying attention to the gauge invariance of the 
probe interaction. 

The internal structure of the resonant A(1405) with full 
coupled-channels has been studied in two ways, one on 
the resonance pole position in the complex energy plane 
and the other on the real energy axis around the A(1405) 
resonance region. The first method has produced the 
resonant A(1405) structure as the matrix clement of the 
probe current by the resonance vector without contam- 
ination from nonresonant background, and exactly kept 
the gauge invariance. The second one, on the other hand, 
has given values which may be observed in experiments 
while the nonresonant contamination has been involved. 
We have found that the resonant A(1405) state has softer 
form factors and larger spatial radii than those of the 
typical baryon, and the structure is largely dominated 
by the KN component. The charge distribution has in- 
dicated that the light K distributes around the proton 
which sits in the central region. The characteristic struc- 
ture of the A(1405) shown in the resonance form factor 
on the resonance pole position can be kept in the effective 
form factor on the real axis. We have also found that the 
decay of the A(1405) into the ttE channel through the 
photon coupling causes an escaping oscillation pattern in 
the density distribution in the coordinate space. 

We have also discussed internal structure of a KN 
bound state without the decay width to extract an in- 
tuitive picture of the bound state. We have found the 
similar structure with the resonant A(1405), thanks to 
the dominance of the KN component. We also have ob- 
served that such a structure shrinks (stretches) as the 
binding energy increases (decreases), which is consistent 
with the behavior of the bound state in the quantum 
mechanics. We have seen the meson mass dependence 
of the bound state which indicates that the behavior of 



the relative distribution of the meson and nucleon in- 
side the bound state is consistent with the kinematics of 
meson-nucleon system. These results have verified that 
the form factor defined through the scattering amplitude 
serves as a natural generalization of the form factor for 
the resonance state. Also the relation between the cou- 
pling strength and the meson mass shows that the physi- 
cal kaon mass appears to be within the suitable range to 
form a molecular bound state with the nucleon through 
the chiral SU(3) interaction. 
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Appendix A: Normalization of form factors 

In this Appendix we discuss the normalization of the 
MB-f* M'B' amplitude Tj!^{^; Q^) induced by a 
conserved current in the soft limit — 0, where 
y/s and are the meson-baryon center-of-mass energy 
and the squared momentum transfer with opposite sign, 
respectively. Because we are considering the matrix el- 
ement with the conserved current J^, the amplitude 
MB^* — > M'B' should automatically be normalized by 
Ward-Takahashi identity: 

Q^^^^^ = T,%"iV:s;Q'^0), (Al) 

where T(Vs) is the MB ^ M'B' amplitude, Q the 
"charge" of the meson-baryon system with respect to the 
probe current, and the indices i and j stand for the final 
and initial meson-baryon channels, respectively^^. Thus, 
we prove that the amplitudes discussed in sections II C 
and II D satisfy Eq. (Al), concentrating on the — Q 
case. 



It is possible to generate a resonance, such as the a meson in 
TTTT scattering, the lower pole of the A(1405) resonance in ttE 
scattering, and the a\ meson in 7rp scattering. 



The Ward identity of the general unitarized amplitude was dis- 
cussed in Ref. [64]. Here we confirm the Ward-Takahashi identity 
in a derivative form (Al) within our formulation. 
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In this study the meson-baryon amphtudc T is ob- 
tained as the solution of the Bethe-Salpeter (BS) equa- 
tion 



7 ^ VikGkTkj — Vij 



J^T.kGkVkj, (A2) 



with the interaction kernel V and the meson-baryon loop 
integral G [see Eq. (33)]. The amplitude is obtained 
by evaluating ten diagrams given in Figs. 3 and 4 in which 
the probe current is attached to the meson-baryon scat- 
tering amplitude. Among them, the amplitudes Tji^^^y 

and T^i^r^^ diverges in the soft limit — > 0. 



7(8)' ^7(9) 



7(10) 



These terms correspond to the infrared divergences com- 
ing from the bremsstrahlung processes of the external 
hadrons and are canceled by the vertex corrections, leav- 
ing no contributions to the amplitude at = 0. 
Hence, in the soft limit, only the six diagrams arc rele- 
vant for T^, 



7(a) 



(A3) 



with each contribution calculated in our formulation as, 



7(1)*J ~ / . ^'ik^Mi.-^kj, 



(A4) 



rpfJ. 

7(2)JJ 



k 



rpfJ. 

7(3)y 



Y^T.kGkVi^iGiTij, 

k.l 



7(4)ij 



J, /J 



rpfJ. 



(A6) 

(A7) 
(A8) 

(A9) 



where the loop integrals with photon couplings D!^, ZJg, 
and the vertex function are defined in Sec. II. 

The key relations for proving Eq. (Al) arc the deriva- 
tive of the vertex and the loop integrals at = 0, 



.dV,j _ I 



, dGk 



(AlO) 



(All) 



Equations (AlO) and (All) are the Ward-Takahashi 
identity for the elementary vertex V and two-body free 



propagator G. Thus, if the vertex T and the loop inte- 
grals Dm and Db conserve the charge associated with the 
current, Eqs. (AlO) and (All) should be satisfied. Then, 
using these relations and the BS equation (A2), we can 
prove Eq. (Al) by an algebraic calculation. Omitting the 
superscript /i = and evaluating all the terms at = 0, 
we can write the left-hand side of Eq. (Al) as, 

[Left-hand side of (Al)] 
dT ^ d 



Q 



d^s 
T + TGT 



Q iy + VGT) 

dy/s 

+ QV 



dG 



T + G 



dT 



dy/s dy/s ) 



(A12) 



where we have use Eqs. (A2) and (AlO). From Eqs. (A7) 
and (A8), the first and second terms in the right-hand 
side of (A12) are T^{4) and T^{5}, respectively. Replacing 
V in the third term of Eq. (A12) with T~TGV, we have. 



[Left-hand side of (Al)] 

= E T^(a) + QiT~TGV) 



a=4,5 



dG 
dyG 



T + G 



dT \ 
d^s) 



a=4,5 



QTG 



- dG 

dy/s 



QTGV^T 



dy/s 



dT 
dy/s 



QTGVG 



dT 
d^/s' 



(A13) 



(A5) From Eqs. (A4) and (A5) with (All), the second term 



of the right-hand side in (A13) is equal to T'^(i) + T^(2)- 
Then, using the product rule of the differentiation for the 
third and fifth terms of Eq. (A13), we have, 



[Left-hand side of (Al)] 



a=l,2,4,5 



[VGT) 



dV 

dy/s 



GT 



QTG 



dT 
dyfs 



(A14) 



From Eqs.(A6) and (AlO), the second term in the paren- 
thesis in Eq. (A14) is T^{^)- Therefore, collecting re- 
mained terms, 

[Left-hand side of (Al)] 

= E T^(a)+QTG-^{T-VGT) (A15) 

Then, using Eqs. (A2), (A9), and (AlO), we finaly achieve 
from the left-hand side of Eq. (Al) to its right-hand side, 

[Left-hand side of (Al)] 

6 

= E^'r(°) " [Right-hand side of (Al)]. 



(A16) 



a=l 



Thus, with (A3), wc obtain Eq. (Al). We note that this 
relation is valid at any value of the meson-baryon centcr- 
of-mass energy y/s, even in complex energy plane. Since 
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Eq. (Al) is the relation between meson-baryon scattering 
amplitude and that with photon attached, the existence 
of the resonance pole is not necessary for the argument 
above. 

The relation (Al) immediately leads to correct normal- 
ization of the effective form factor defined in Eq. (21), as. 



^cff^=0(Q2 ^ 0; ^^) = 



dTij/dy/s 



Q. (A17) 



^=0 



This normalization is achieved at any energies ^/s, re- 
gardless of existence of the resonance, as long as one 
includes all the relevant contributions in Eq. (A3). If 
there exists the resonance pole at the complex energy Zn 
in the meson-baryon scattering, taking limit ^/s — > 
we see that F"^ = Ti^=° / {dT^j / d^/s) coincides with the 
resonance form factor F evaluated in Eq. (19). Hence 
the relation (Al) also guarantees the correct normaliza- 
tion of the resonance form factor F. In the case that 
the energy ^/s is close to the resonance pole position 
{'/s ~ 1420 MeV in our case) on the real axis, the charge 
Q will be dominated by the contributions from the three 
diagrams in Fig. 3, as discussed in Sec. II E. 

We also derive, from Ward-Takahashi identity (Al), a 
relation between the coupling strengths of the resonance 
state to the meson-baryon channel g and the derivative 
of the vertex V and the loop function G: 



El dGi 



dy/S 



-1. (A18) 



To prove Eq. (A18), we consider explicit forms of Tj 
and dTij/dy/s for the resonance contribution. The 



derivative of the MB — > M'B' amplitude dTij/d^/s is 
shown in Eq. (22). The MB7* ^ M'B' amplitude r^=° 
has double-pole terms for the resonance as well as the 
less singular terms. The former corresponds to Fig. 3 and 
is expressed in Eqs. (A4)-(A6), while the latter does to 
Fig. 4 and is in Eqs. (A7)-(A9). Using the explicit form 
of Tij given in Eq. (11) and the Ward-Takahashi iden- 
tity (AlO) and (All) for the double-pole terms, T^=° 
can be expressed as. 



0) 
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k.l 



gi9k 



dCk 



( 



>kl 



r< dVki 

dy/S 



9i9j 



less 



(A19) 



where T^^^*^ represents the less singular terms than 

the double-pole contribution. Thus, calculating Q = 
T^=° / {dT^jl dy/s) and taking hmit Zr, both the 

nonresonant contribution in dTij / d^ and less singular 
terms T}^^-^ in T^~^^ automatically drop and we have, 



Q = 



T: 



dTij/dy/s 



gugi 



k,l 



\d^/s 



dy/s 



(A20) 



hence, we obtain Eq. (A18). This is the exact form of 
the relation between the coupling strength gi and the 
derivative of the loop integral dCi/d^/s, discussed in [30, 
69]. Equation (AlS) exactly corresponds to the relation 
known as the Ward identity: 



^(1)^(2) 



(A21) 



with the vertex-renormalization factor Z^^^^ and the wave- 
function one Z(^2)i respectively. In our case with the 
resonance state, Z(2) is the residue of the resonance 
contribution in the meson-baryon scattering amplitude, 
Z{2)ij ~ ~9igj7 E^nd Z,~J is the current-coupling term. 



^(1) 
fdC,^ 



d^/s 



(A22) 



through Eqs. (AlO) and (All). Using -2^(1)^ and Z(2)ij, 
Eq. (AlS) can be rewritten as, 



tr 



^(1) ^(2) 



(A23) 



Appendix B: Loop integrals and 

In this Appendix we discuss the properties of the loop 
integral -Df^ 3(^/3; Q^) [given in Eqs. (52) and (53)], 
which describes that one photon is attached to the one of 
the internal propagators. For concreteness, we consider 
the electromagnetic probe interaction for the A(1405). 

Integrating over the four-momentum qf^ in Eqs. (52) 
and (53) in the Breit frame for the A(1405), these func- 
tions are explicitly written with the Feynman parameters 
X and y as 
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FIG. 23: The real part (left) and the imaginary (right) parts of the loop integrals D^j + of the K p channel (upper) and 
the 7r"'"E~ channel (lower) as functions of \m[z\ and . The real part of the energy is chosen to be Re[2;] = 1420 MeV. 
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2(,s + y(l-y)Q2) 



+ 



4g2s _ 4j/(l - y)Q'ml - (s + y(l - j/)Q2)(, 



2(s + y(l - y)02)^4g2s - 4y(l - is)Q'^rai + ze 



-L (v^; Q^; m^, A4) 



(B2) 
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(47r)2 



dx / dy 



(47r) JO JO 
fik f icr X q 



2Af„ 



In 



xMl + (1 — x)m\ — x{l — x)s + — y)Q'^ 



icr X q 

2M„ 



dy- 



2Ml 



y{i-y)Q^ 



Ml + y[l-y)Q^ 



(s + ml 



Ml) 



- 42/(1 - y)Q'^ml + i 



:L Q2. jn^^ Mk) 



(B3) 



L (Vi; 0'; A4, mk) ^In .s + il/^ - ml + J Aqls - 4y(l - y)QHll + ze 



In ,s - Mi + + 2y(l - zj)g^ + . /4g2,s _ Ay{l - y)Qni, 



-In -s-Ml 



Aqls - 4y(l - y)Q^Ml + 



In 



+ Ml - - 2y(l - y)Q2 + ^4g2s - 4y(l - y)Q2M2 + ze 
I 



(B4) 



Here Qm^i Qb^^ ^-nd /^b^ are the charge of the meson 
and baryon, and magnetic moment of the baryon in k 
channel, respectively, and qk is defined in Eq. (29) so 



that 4q|s is simply written as. 



A(s, m 



I Ml) 



(B5) 



with the Kiillen function A(x, y, z) — x +y + z —2xy — 
2yz — 2zx. We note that spatial components of the loop 
integrals in which photon couples to mesons vanish in the 
Breit frame, D'^ = (a = 1, 2, 3), since the pseudoscalar 
mesons do not have the magnetic moments. In the above 
equations the x integrations are analytically performed; 
the y integrations can in principle be done analytically, 
but we perform the numerical integration in the practical 
calculation with keeping ie to be a small finite value, 
which ensures the correct boundary condition during the 
analytic continuation. 

As one can see from Eqs. (B1)-(B4), the integrands 
in g are obtained in the analytic form for the energy 
^/s and the squared photon momentum Q^, hence we can 
make an analytic continuation from the real to the 
complex value, y/s — ^ z, for the integrands. 

In our study we perform the analytic continuation of 
g in a following way. At first, we fix = and make 
y/s a complex value, -^/s — > y/s + ib = z (6: real) so that 
B continuous with respect to b. In this condition, 
we are in a first (second) Riemann sheet with respect to 
with 6 > (6 < 0) if the energy y/s is above the 
meson-baryon threshold. Then, we make vary from 
to finite real values so that D'^ g are continuous with 
respect to Q^. 

In order to visualize this, we plot the sum of the 
loop integrals + ZJg, which appears in the ampli- 



and TT+E^ channels with the real part of the energy be- 
ing fixed as Re[z] = 1420 MeV, around the pole posi- 
tion of the A(1405). Here we note that in all meson- 
baryon channels for total charge-zero system, the func- 
tions -I- vanish at = at any complex energies 
z, because at = the relations (54) and (55) lead to 

-I- 



iEM 



0. 



Now we discuss the behavior of the K^p loop integral. 
One notice that Im[D^ -I- Dg] = with real-valued en- 
ergy, that is, Im[z] = 0. This is interpreted as the fact 
that K~ and p cannot go on-shell simultaneously with 
the energy below the K~p threshold, y/s < m^- + M-p. 
With the real-valued energy (Im[2;] = 0), the real part 
Re [13^ + ZJg] decreases as we increase . This prop- 
erty of decreasing Re[Z)^ -I- Dg] is visible until around 
Im[z] ~ ±20 MeV in the analytic continuation to imag- 
inary energy, ^/s — >■ z. For the K^p channel, D^^ + _Dg 
is smooth in the Q^-Im[z] plane. 

The 7r+S~ channel of the loop integral, on the other 
hand, has some nontrivial behaviors. One is the diver- 
gence of the loop integral at finite with real ^s, which 
is clearly seen in the imaginary part ^^[Z?^ -I- Dg], and 
another is the discontinuity in the real part Re[£'^ -I- ZJg] 
running above the same . The divergences of Dy^ and 
Dg take place, as will be discussed later, at 



Ml 



(B6) 



respectively. Hence, the loop integral with > di- 
verges only when the real-valued energy ^/s is above the 
two-body threshold, so that Aq^s > 0. Furthermore, this 
divergence point corresponds to the t-channel threshold; 
with fixed -y/i above the threshold, all the three prop- 
agating particles can simultaneously go on-shell in Dm 
tude as in Eqs. (48)-(50), in the Q^-Im[z] plane. (Dg) for < Aq'^s/M'^ (Q^ < Aq'^s/m'^), while for 



In Fig. 23 we show time component D' 



D'^ in K-p 



> Aq^s/M^ [Q^ > Aq^s/m^) two propagating mesons 
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(baryons) cannot simultaneously go on-shell in Dm (Db) 
if we constrain the baryon (meson) to go on-shell. 

The origin of the divergences of Dm, and Db^. is the 
factor 1 / y/4qf.s — 4y(l — y)Q^fnj. + ie in the integrands 
of Eqs. (B1)-(B3), where fhk represents the meson mass 
TOfc (the baryon mass Mjt) when the photon couples to 
the baryon (meson) propagator. By dropping the ir- 
relevant factors, which do not cancel singularity from 
1 / ■\/4g^s — 4?/(l — y)Q^m| + ie, and performing the y 
integration, we obtain. 



respectively. Therefore, I^^Ztxt. ^^'^ limit. 



It ^ I dy 



2Qmk 



In 



V4g2s - 4j/(l - y)Q^ml + ie 



(B7) 



which contains the logarithmic divergence at Qfhk = 
±2qk\/s- Therefore, the loop integral Dm^, {DBk) log- 
arithmically diverges at = Aq^s/Mk {Q^ = Aq\s/mk) 
with ^/s > TOfe -|- Mk ■ 

Now it is instructive to investigate the behavior of /^'^ 
with respect to Im[2] and Q^. In Fig. 24 we plot I^ttS 
a function of Im[z] and with conditions fhk = M-^ and 
Re[z] = 1420 MeV. Here analytic continuation is per- 
formed in the same way as for Dm and Db . From Fig. 24 
we can see the real part of Ik^TrS contains divergence 
point at Im[z] = and = Aql^s/Ml ~ 0.16 GcV^ 
whereas the imaginary part of itiri: shows discontinu- 
ity along the real energy axis above the same value 
(~ 0.16 GeV^). These properties are explained by us- 
ing analytic properties of a logarithm function ln(a;). 
Namely, taking limit (2q'7rS y/s — QMy) from any 
complex energy z and any real Q, we have, 

ln(2g^s - QM^) = In \2q^^ ^/i - QM^\ + iO, (B8) 

with argument of (2g7rS \/s — QM^), 9, which is finite, 
hence the divergence appears only in the real part of 
-^fc=7rS- Further, in order to see the behavior of the 
imaginary part of itirY.^ ^^^^ energy with 

an infinitesimal imaginary part as ^/s ± ie (e > 0). 
Then, from y/s the divergence point of itni: fixed as 

1/2 

Qo = (v^)/-^^ with the square root of the Kallen 
function ^J-^{y/s) — 2q-^^y/s. Now let us consider Ik=Tri: 
with the energy ^/s ± ie and momentum Q — Qo + AQ, 
where AQ is finite positive value. Since the Kallen func- 
tion gives a relation. 



\'J^iV~s±ie)^X'J^{^s)±iS, 



(B9) 



where S is an infinitesimal positive value determined by 
the \/s and e and has a property lime_j.o S ~ 0, the numer- 
ator and denominator inside the logarithm in Eq. (B7) 
become: 

X'J^iV^ ± ie) + QMs = (2Qo + AQ)A/s ± iS, (BIG) 
ii/^r /o fl r„ — Ani\j^j-nX (Bll) 



•In 



/(2Qn + AQ)A/s±i(5 



2Q, 



-AQAfs ± iS 
+ AQ 



AQ 



=F iTT (e, S 0), 



(B12) 



hence the discontinuity appears in the imaginary part of 

7"div 

1 I „v ■ 



1 



As a consequence of the factor in the integrands, 
42/(1 — y)Q^'fnk + ie, Dm, and Db, in the 



Xj^iV^ ± ie) ~ QMs = -AQMs ± iS, 



channels with energy above the threshold contains di- 
vergence in the imaginary part at Q^ = Aq'^s/mf. and 
the discontinuity in the real part along the real energy 
axis above the same Q^ value. Here we note that the L 
function (B4) is pure imaginary so that the discontinuity 
(divergence) of the imaginary (real) part of /^'^ is man- 
ifested in the real (imaginary) part of the loop integrals 

DMfc,Bfc- 

These diverging structures in the loop functions Dm,b 
eventually lead to the divergences of the effective form 
factors (21) at Q^ = Aifs/fn^, and to the peak struc- 
tures of the form factors (19) around that Q^. Then, the 
"peaks" of the (effective) form factors in the momentum 
space lead to the oscillation components in the density 
distributions in the coordinate space with period given by 
2-K / {2qs/s / rfi) = ■Kfh/{qy/s) through the Fourier trans- 
formation. Therefore, for the channel whose threshold 
is lower than the total energy y/s, we will have an os- 
cillation component of the "decaying part" , which can 
be interpreted as the decay of the system into the open 
channels with kicked meson and baryon to the on-shell 
through the photon coupling. 

Here let us show the behaviors of K~p and 7r+S~ com- 
ponents of effective form factor and density distribution 
in electric probe with energy fixed as ^/s = 1420 MeV. 
The form factor and density distribution on the reso- 
nance pole position are shown in Sec. HI A. As you 
can see, the effective electric form factor in tt+S^ chan- 
nel shows divergence at Q^ = Aq^^s/M^ ~ 0.16 GeV^, 
which generates an oscillation behavior in the effective 
electric density distribution with the period close to 
^Ms/(gyi) ~ 3.1 fm. 

Finally we note that the oscillation part will be clearly 
seen only in the case that the probe current couples to 
the mesons. In this case m corresponds to the baryon 
mass M, hence the value Aq^s/JVP appears around typ- 
ical hadronic scale Q^ < IGeV^, which will be impor- 
tant for the dynamics of the hadrons. This is not the 
case when the probe current couples to baryons, since in 
this case the value Acp's/m? becomes much larger than 
1 GeV^, due to the light meson masses to, especially for 
pions, 771 ~ m^r, and such a high Q^ coupling should 
be strongly suppressed by form factors of constituent 
hadrons. 
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Re. part 




Im. part 




FIG. 24: The real part (left) and the imaginary (right) parts of the /fei^^s (B7) as functions of Im[2;] and . The real part of 
the energy and fhh are chosen to be Re[z] = 1420 MeV and fhh = Ms, respectively. 




FIG. 25: Meson-baryon components of the effective electric form factors {F^ , left) and effective electric density distribution 
(P|?, right) on the real energy axis. Calculations are performed with the center-of-mass energy y/s = 1420 MeV. 
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